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The  overarching  goal  of  this  grant  award  is  to  characterize  differential  splicing  of  oncongenes  and  tumor  suppressor  genes  in  prostate 
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signaling  protein  as  defined  by  in  vitro  assays  and  mouse  xenograft  studies,  ii)  is  associated  with  worse  prognosis  in  patients,  and  iii)  is 
resistant  to  small  molecule  inhibitor  idelalisib  or  CAL-101  (selective  PI3KC8  inhibitor;  FDA  approved  for  treating  hematologic  cancers). 
Our  discovery  portends  a  genetic  screening  test  for  patients  with  aggressive  tumors  that  are  resistant  to  idelalisib  therapy.  In  year  3  of  this 
award,  we  have  focused  our  efforts  on  characterizing  the  molecular  mechanism  underlying  increased  oncogenic  signaling  by  PI3KC8-S 
compared  to  PI3KC8-L  (expressed  in  CA  PCa).  Our  findings  demonstrate  that  PI3KC8-L  preferentially  complexed  with  the  regulatory 
subunit  p85a,  leading  to  decreased  kinase  activity  compared  to  PI3KC8-S  that  has  lower  affinity  for  p85a.  PI3KC6-L,  in  the  absence  of 
p85a,  exhibited  kinase  activity  comparable  to  PI3KC8-S  but  was  still  sensitive  to  the  inhibitory  effects  of  idelalisib. 
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1.  INTRODUCTION 

There  are  striking  population/race  disparities  in  prostate  cancer  (PCa)  risk  and  survival  outcome 
borne  out  of  current  health  statistics  data.  This  is  particularly  evident  between  African  Americans 
(AA)  and  their  European  American  (EA)  counterparts.  Epidemiologic  studies  have  shown  that 
higher  mortality  and  recurrence  rates  for  prostate  cancer  are  still  evident  in  AA  men  even  after 
adjustment  for  socioeconomic  status,  environmental  factors  and  health  care  access.  Thus,  it  is 
likely  that  intrinsic  biological  differences  account  for  some  of  the  cancer  disparities.  Our 
overarching  hypothesis  is  that  the  biological  component  of  prostate  cancer  health  disparities  is 
due,  in  part,  to  population-dependent  differential  splicing  of  oncogenes  and  tumor  suppressor 
genes  in  cancer  specimens.  The  application  of  genomic  approaches  has  identified  splice  variants 
in  AA  specimens,  but  absent  in  EA  specimens,  encoding  more  aggressive  oncogenic  proteins, 
thereby  producing  a  more  cancerous  phenotype. 


2.  KEYWORDS 

Prostate  cancer,  cancer  health  disparities,  alternative  RNA  splicing,  African  American,  European 
American,  oncogenes,  tumor  suppressor  genes,  phosphatidylinositol-4,5-bisphosphate  3-kinase 
catalytic  subunit  delta,  fibroblast  growth  factor  receptor  3 


3.  ACCOMPLISHMENTS 
Year  3  goals  as  stated  in  SOW: 

Specific  Aim  1.  To  define  splice  variant  pairs  (AA-specific  variant  versus  EA-counterpart 
variant)  associated  with  differential  oncogenic  behavior  in  vitro,  and  to  delineate  the 
mechanism  of  action. 

Task  1.  Full-length  cloning  and  in  vitro  validation  of  splice  variant  pairs.  Subtasks  will  be 
run  concurrently  and  are  as  follows: 

la.  Full-length  cloning  of  splice  variant  pairs  and  ectopic  over-expression  into  PCa  cell  lines. 

lb.  In  vitro  validation  of  differential  oncogenic  behavior  by  full-length  splice  variant  pairs. 
Splice  variant  pairs  (e.g.  AA-specific  versus  EA-counterpart  variant  of  PIK3CD  and 
FGFR3)  will  be  individually  over-expressed  in  the  same  PCa  cell  line  background,  and 
screened  for  differential  oncogenic  behavior. 

lc.  In  vitro  validation  of  differential  protein/enzyme  activity  by  full-length  splice  variant 
pairs.  Splice  variant  pairs  will  be  individually  over-expressed  into  appropriate  cell  line  for 
enzyme  activity  assays  and/or  assessment  of  downstream  activation  of  cell  signaling 
components.  Activation  of  downstream  signaling  components  by  splice  variants  will  be 
assessed,  for  example,  by  measuring  phosphorylation  of  downstream  signaling 
components  with  phospho-specific  antibodies  (e.g.  phospho-Akt,  phospho-ERK,  etc.). 

Task  2.  In  vitro  screening  and  full-length  cloning  of  additional  splice  variant  pairs.  Sub  tasks 
will  be  run  concurrently  and  are  as  follows: 

2a.  Exon-targeting  and  splice  junction- targeting  siRNAs  will  be  used  in  appropriate  PCa  cell 
lines  to  identify  splice  variant  pairs  exhibiting  differential  oncogenic  behavior  following 
knockdown. 
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2b.  From  subtask  2a,  we  will  select  5-10  splice  variant  pairs  that  exhibited  differential 

oncogenic  behavior  for  full-length  cloning  and  ectopic  over-expression  in  appropriate  cell 
lines. 

2c.  Cell  lines  over-expressing  individual  full-length  variant  pairs  (e.g.  AA-specific  variant 
versus  EA-counterpart  variant)  will  be  validated  in  vitro  for  differential  oncogenic 
behavior  using  in  vitro  screens  described  in  sub  task  lb.  We  will  also  test  for  differential 
sensitivity  of  splice  variant  pairs  to  small  molecule  inhibitors,  if  available. 

2d.  Cell  lines  over-expressing  individual  variant  pairs  (e.g.  AA-specific  variant  versus  EA- 
counterpart  variant)  will  be  screened  in  vitro  for  differential  protein/enzyme  activity  and 
cell  signaling  as  described  in  sub  task  Ic.  We  will  also  test  for  differential  sensitivity  of 
splice  variant  pairs  to  small  molecule  inhibitors,  if  available. 

Specific  Aim  2,  To  characterize  oncogenic  differences  of  splice  variant  pairs  in  vivo  using 
xenograft  animal  models. 

Task  1.  Validate  differential  oncogenic  behavior  of  the  splice  variant  pair  for  PIK3CD  in 

vivo.  Stably  expressed  S  (AA-specific)  or  L  variants  (EA-counterpart)  of  PIK3CD  in  appropriate 
cell  line(s)  will  be  transplanted  (1x10^  to  10^  cells)  into  male  SCID-NOD  immuno-deficient  mice 
for  proliferation  and  metastasis  assays. 

Task  2.  Validate  differential  oncogenic  behavior  of  additional  splice  variant  pairings  in  vivo. 

We  will  test  in  vivo  an  additional  4-9  splice  variant  pairings  defined  in  Aim  1,  Task  1,  Subtasks 
Ib-lc  (e.g.  one  variant  pairing  could  be  the  AA-specific  and  EA-counterpart  variants  for 
FGFR3),  or  defined  in  Aim  1,  Task  2,  Subtasks  2c-2d. 


Year  3  major  accomplishments  include  the  following; 

i.  Manuscript  detailing  our  findings  on  the  oncogenic  behavior  of  the  AA- 
specific/enriched  PIK3CD  short  variant  (PIK3CD-S)  compared  to  the  EA  PIK3CD 
long  variant  (PIK3CD-L)  has  been  positively  reviewed  by  Nature  Communications 

Assessment  of  our  Nature  Communications  manuscript  by  the  reviewers  was  positive,  albeit  they 
requested  a  number  of  additional  experiments  that  have  been  completed  this  past  month.  Our 
manuscript  has  been  resubmitted.  Our  additional  findings  can  be  summarized  as  follows: 

AA-specific/enriched  variant  PI3KCD-S  is  resistant  to  small  molecule  inhibitor  CAL-101, 
while  EA  variant  PI3KCD-L  is  sensitive.  In  our  original  Nature  Communications  submission, 
we  had  stably  over-expressed  the  PIK3CD-S  and  PIK3CD-L  variants  (individually  in  the  EA 
PCa  cell  lines  PC-3  and  VCaP.  These  stably  transfected  cell  lines  were  tested  in  vitro  for 
sensitivity  to  CAL-101  treatment  (Fig.  lA).  CAL-101  is  a  PI3K5  inhibitor  in  clinical  trials  for 
various  cancers.  PCa  cell  lines  over-expressing  the  EA  P1K3CD-L  variant  exhibited  a  decrease  in 
the  activity  of  the  PI3K/AKT  pathway  following  CAL-101  treatment,  as  seen  by  a  loss  of  AKT, 
mTOR  and  S6  phosphorylation.  Remarkably,  the  same  EA  PCa  cell  lines  stably  over-expressing 
equivalent  levels  of  the  AA  P1K3CD-S  variant  were  completely  resistant  to  CAL-101 .  In  other 
words,  there  was  no  significant  change  in  AKT,  mTOR  and  S6  phosphorylation  levels  before  and 
after  CAL-101  treatment. 
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In  cell  proliferation  assays,  we  demonstrate  that  BrdU  labeling  in  VCaP  and  PC-3  cells  over¬ 
expressing  the  EA 
PIK3CD-L  variant  was 
inhibited  by  CAL- 101  in  a 
dose-dependent  manner, 
while  proliferation  of  VCaP 
and  PC-3  cells  over¬ 
expressing  the  AA 
PIK3CD-L  variant  were 
resistant  to  CAL- 101  (Fig. 

IB). 

The  reviewers  had 
requested  that  we  assess 
whether  or  not  cell  lines 
over-expressing  the 
PIK3CD-S  or  PIK3CD-L 
variant  exhibited 
differential  sensitivity  to 
inhibitors  of  downstream 
signaling  components. 

Interestingly,  PIK3CD-S- 
and  PIK3CD-L-OYer- 
expressing  cell  lines  were 
equally  to  selective 
inhibition  of  AKT  with 
MK-2206  (Fig.  1C).  These 
results  demonstrate  that 
PI3K5-S-stimulated 
proliferation  is  resistant  to 
CAL-101  inhibition  in 
sharp  contrast  to  PI3K5-L; 
while  inhibition  of  AKT, 
which  is  downstream  of 
PI3K5-S,  effectively 
blocked  proliferation. 

ii.  Defining 

mechanism  of  PI3KC6-S  oncogenecity  and  resistance  to  CAL-101  inhibition 

Employing  cell-free  system  to  define:  a)  mechanism  of  increased  oncogenecity  of  PI3KC8-S 
over  PI3KC6-L,  and  b)  the  mechanism  of  PI3KC6-S  resistance  versus  PI3KC6-L 
sensitivity  to  CAL-101  inhibition.  A  major  question  that  the  reviewers  of  Nature 
Communications  required  us  to  address  was  the  mechanism  of  enhanced  PI3KC8-S 
oncogenecity  and  resistance  (Fig.  2).  The  result  of  the  exclusion  of  exon  20  (168  bp)  in  the 
P1K3CD-S  variant  is  an  in-frame  deletion  of  56  amino  acids  (residues  810-865)  in  the  catalytic 
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Fig.1  .  PI3K6-S,  but  not  PI3K6-L,  is  resistant  to  smaii  moiecuie  inhibition  of 
PI3K/AKT/mTOR  signaiing  and  proliferation,  (a)  Assessment  of 
PI3K/AKT/mTOR  signaling  following  treatment  with  vehicle  (saline)  or  CAL-101 
(100  nM,  24  hrs).  PI3K/AKT/mTOR  signaling  was  assessed  by  western  blot 
analysis  with  phospho-antibodies  to  AKT  (pAKT)  and  S6  ribosomal  protein  (pS6). 
P-actin  served  as  a  loading  control.  His-tag  antibody  was  used  to  demonstrate  equal 
expression  of  His-tagged  variant  PI3K5  protein  in  stably  transfected  cell  lines. 
Representative  images  from  n  =  3-4  independent  western  blot  experiments. 
*Significantly  different  from  corresponding  vehicle  control,  P  <  0.05  by  ANOVA 
with  Dunnett’s  post-hoc  test,  (b)  Proliferation  in  VCaP  and  PC-3  cells  stably 
overexpressing  the  P1K3CD-S  variant  or  P1K3CD-L  variant  following  treatment 
with  vehicle  (saline)  or  selective  PI3K5  small  molecule  inhibitor  CAL-101  (24 
hrs),  and  (c)  treatment  with  vehicle  (saline)  or  selective  AKT  small  molecule 
inhibitor  MK-2206  (24  hrs).  Proliferation  was  assessed  using  a  BrdU  labeling 
assay.  Data  presented  as  the  mean  ±  SEM  of  n=4-6  independent  experiments  for 
each  treatment  group.  *  or  #Significantly  different  from  corresponding  vehicle 
control,  P  <  0.05  by  ANOVA  with  Dunnett’s  post-hoc  test. 
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domain  of  the  PI3K6-S  isoform  (Fig.  2a).  To  gain  further  insight  into  the  functional  differences 
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Fig.  2.  Cell-free  kinase  assay  of 
PI3K6  isoforms  and  tertiary 
strueture  modeling  of  PIK38- 
small  moleeule  inhibitor 
interaetions.  (a)  Schematic 
representation  of  protein  domains 
of  PI3K5-L  and  -S  isofonns. 
Adaptor  (p85)-binding,  RAS 
binding,  C2,  helix  and  catalytic 
domains  are  highlighted. 
Phosphorylation  and  ubiquitination 
sites  (S,  T,  Y  and  K)  and  the 
region  encoded  by  exon  20  (56 
amino  acids)  residing  in  the 
catalytic  domain  are  indicated,  (b) 
Co-immunoprecipitation  (Co-IP) 
of  PI3K5/p85  complex  followed  by 
western  blot,  (c)  and  PI3K  activity 
assays.  S,  supernatant;  W,  wash 
fraction;  E,  eluted  fraction.  Blots 
are  representative  of  3-5  similar 
results.  Anti-His  antibody  was 
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used  in  the  Co-IP  experiments  and  anti-His,  anti-p85a  and  anti-actin  antibodies  were  used  in  the  western  blotting. 
*Significantly  different  kinase  activities  in  total  lysates  of  PI3K5-S  vs.  PI3K5-L  expressing  cells,  or  purified  PI3K6-S/p85  vs. 
PI3K5-L/p85  complexes.  P  <  0.05  using  Student  t-test.  Data  presented  as  mean  ±  SEM  of  n  =  4  for  each  treatment  group,  (d) 
Purification  of  His-tagged  PI3K5-L  and  -S  isofonns.  Western  blot  analysis  of  Ni-NTA  resin-purified  PI3K5  isoforms  from 
transfected  PC-3  cells  using  His  and  PI3K6  antibodies.  FT,  flow-through;  W,  wash  fraction;  E,  eluted  fraction.  Blots  are 
representative  of  4  independent  experiments  with  similar  results.  Closed  arrowheads  indicate  PI3K5  isoforms,  (e)  Cell-free 
kinase  assay  of  L  and  S  isoforms  of  PI3K6  in  the  presence  of  vehicle  (PBS),  lOOnM  wortmannin  or  lOOnM  CAL-101. 
*Significantly  different  from  vehicle  control-treated  PI3K6-L  isoform.  P  <  0.05  by  ANOVA  with  Dunnett’s  post-hoc  test. 
Data  presented  as  mean  ±  SEM  of  n  =  4-5  for  each  treatment  group. 


between  PI3K8  isoforms,  the  interaction 
of  PI3K8-L  and  -S  with  regulatory  subunit 
p85a  was  investigated.  Whole  cell  lysates 
from  PC-3  cells  over-expressing  either 
His-tagged  PI3K8-S  or  PI3K8-L  were 
subjected  to  western  analysis, 
demonstrating  that  each  cell  line  expressed 
equivalent  levels  of  their  respective  PI3K8 
isoform  as  well  as  equal  p85a  expression 
(Fig.  2b;  left  panel).  Interestingly,  co-IP  of 
the  PI3K8/p85a  complex  from  whole  cell 
lysates  using  an  anti-His  antibody 
demonstrated  that  p85a  bound  with  3-4- 
fold  greater  proficiency  to  PI3K8-L 
compared  to  p85a  binding  to  PI3K8-S 
(Fig.  2b;  right  panel,  column  E).  Binding 
proficiency  was  inversely  correlated  with 
PI3K8  isoform  kinase  activity  (Fig.  2c; 
right  panel). 


Fig.  3.  Knockdown  of  p85a  increases  invasive  activity  of  EA 
PCa  cell  lines.  Invasion  was  assessed  using  Matrigel  assay 
following  siRNA-mediated  knockdown  of  p85a  (PIK3R1). 
SiRNA-SMART  pool  for  targeting  PIK3R1  was  purchased  from 
Dhannacon  (Lafayette,  CO).  Data  presented  as  the  mean  ±  SEM 
of  n=3-4  independent  experiments  for  each  treatment  group.  *P 
<  0.05  by  ANOVA  with  Dunnett’s  post-hoc  test.  siNS, 

Nonsense  siRNA  control. 
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Next,  His-tagged  PI3K6  isoforms  were  purified  from  the  lysates  of  PC-3  cells  overexpressing 
either  PI3K6-S  or  PI3K6-L  using  Ni-NTA  resin  columns.  As  shown  in  Fig.  2d  (left  and  middle 
panels),  purification  of  the  PI3K8-S  and  -L  isoforms  was  verified  by  western  blotting  using  anti- 
His  or  anti-PI3K8  antibodies.  In  addition,  the  Ni-NTA  resin  column  purification  approach 
resulted  in  the  isolation  of  PI3K8  isoforms  that  were  no  longer  bound  to  p85a  (Fig.  2d;  far  right 
panel,  column  E).  The  purified  PI3K8  isoforms  (minus  p85a)  were  incubated  with  vehicle,  non- 
selective  PI3K  inhibitor  wortmannin  (100  nM)  or  PI3K8-specific  inhibitor  CAL-101  (100  nM), 
and  subjected  to  a  PI3K  activity  assay.  In  the  absence  of  bound  p85a,  kinase  activity  of  PI3K8-L 
was  equivalent  to  PI3K8-S  (Fig.  2e,  compare  vehicle  treatments).  In  agreement,  siRNA- 
mediated  knockdown  of  p85a  in  wild-type  EA  PCa  cell  lines  VCaP  and  PC-3  was  associated 
with  an  increase  in  invasive  activity  (Fig.  3).  Remarkably,  wortmannin  and  CAE-101 
significantly  inhibited  the  activity  of  the  PI3K8-E  isoform,  but  not  the  PI3K8-S  isoform  (Fig. 
2e).  These  results  demonstrate  that  PI3K8-S  maintains  kinase  activity  even  in  the  presence  of 
small  molecule  inhibitors,  supporting  our 
earlier  in-vitro  and  in-vivo  results 
presented  in  our  Year  2  progress  report. 

iii.  We  have  determined  that  high 
PI3KCD-S  and  FGFR3-S 
expression  is  associated  with 
poorer  patient  prognosis 

We  have  performed  survival  analysis  of 
The  Cancer  Genome  Atlas  (TCGA)  RNA- 
Seq  data,  demonstrating  that  either  high 
expression  of  PIK3CD-S  or  FGFR3-S 
associates  with  significantly  poorer 
survival  in  PCa  patients  (Fig.  4), 

demonstrating  further  the  potential  relevance  of  these  two  short  variants  in  PCa. 


Fig.  4.  PIK3CD-S  (A)  and  FGFR3-S  (B)  are  associated 
with  decreased  survivai  among  PCa  patients.  RNA- 
Seq  and  disease  free  survival  data  for  PCa  (n  =  494 
patients)  was  obtained  from  The  Cancer  Genome  Atlas 
(TCGA)  (https://tcga-data.nci.nih.gov).  P-values  for 
survival  curves  determined  by  the  log-rank  test. 


Year  3  opportunities  for  training  and  professional  development: 

AACR  Health  Disparities  Presentation.  Year  3  of  this  proposal  has  provided  hands-on  training 
for  PhD  graduate  student  Jacqueline  Olender.  The  PI  is  serving  as  Ms.  Olender’s  mentor  and  she 
has  participated  in  both  the  in  vitro  and  in  vivo  work  described  herein.  This  work  is  part  of  Ms. 
Olender’s  PhD  dissertation  research  project.  An  NIH/AACR  Scholar-in-Training  Award  was 
awarded  to  Ms.  Olender  to  present  her  thesis  work  entitled,  “Alternative  splicing  of  FGFR3  as  a 
mechanism  for  prostate  cancer  health  disparities”  at  the  9th  AACR  Conference  on  The  Science 
of  Cancer  Health  Disparities  in  Racial/Ethnic  Minorities  and  the  Medically  Underserved. 


Dissemination  of  results  and  outreach  to  communities  of  interest: 

Mentoring.  During  Year  3  of  this  grant,  my  lab  has  mentored  4  medical  students  (Eauren 
Simpkins,  Daniel  Szvarca,  Breeana  Johnq,  Erin  Good)  during  the  summer  of  2016,  on  various 
PCa  disparities -related  research  projects,  under  the  supervision  of  myself  and  Dr.  Wang. 
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Invited  Presentation.  This  work  was  presented  at  the  June  2016,  Khon  Kaen  University  Cancer 
Summit  Meeting  in  Khon  Kaen,  Thailand.  “Aberrant  Alternative  Splicing  of  mRNA  as  a 
Mechanism  in  Prostate  Cancer  Disparities”,  Lee  N.H. 


Goals  for  next  reporting  neriod: 

None  to  report 


4.  IMPACT 

Impact  on  the  development  of  the  principal  discinlinefsl  of  the  project; 

Principal  discipline  —  Understanding  prostate  cancer  biology  and  disparities.  Taken  together,  our 
in  vitro  and  in  vivo  research  on  A  A- specific/enriched  PIK3CD-S  and  FGFR3-S  variants  provide 
strong  evidence  that  differential  splicing  may  play  a  critical  role  in  PCa  health  disparities.  Our 
future  goal  is  to  identify  the  underlying  molecular  mechanism(s)  responsible  for  the  differential 
splicing  events  observed  in  AA  versus  EA  PCa  specimens. 

Impact  on  other  disciplines: 

Other  disciplines  —  Our  work  demonstrates  that  the  short  variant  encodes  a  more  aggressive 
oncogenic  signaling  protein  isoform,  PI3K8-S,  that  is  resistant  to  small  molecule  inhibitor 
idelalisib  as  defined  by  in  vitro  assays  and  mouse  xenograft  models.  In  contrast,  the 
corresponding  EA  isoform  (PI3K6-E)  encodes  a  less  aggressive  isoform  that  is  sensitive  to 
idelalisib  inhibition.  Our  discovery  portends  a  genetic  screening  test  for  cancer  patients  that  may 
be  resistant  to  idelalisib  therapy,  which  is  particularly  germane  in  hematologic  malignancies. 
Idelalisib  is  EDA  approved  for  relapsed  chronic  lymphocytic  leukemia  and  relapsed  non- 
Hodgkin  lymphoma.  Treatment  resistance  is  observed  in  ~30-40%  of  patients,  and  the 
mechanism  of  resistance  is  currently  unknown.  Our  hypothesis  is  that  the  short  mRNA  variant  of 
PIK3CD  may  be  playing  a  major  role  in  treatment  resistance.  Clinical  and  molecular  studies  are 
currently  underway  to  test  this  hypothesis. 

Impact  on  technology  transfer: 

Our  findings  that  the  AA  PIK3CD-S  variant  protein  is  resistant  to  CAE- 101  has  sparked  interests 
in  companies  that  are  investigating  small  molecule  inhibitors  of  kinases  involved  in  cancer 
progression.  These  companies  are  gaining  an  appreciation  that  alternative  splicing  in  kinases  can 
affect  the  sensitivity  these  signaling  proteins  to  cancer  therapeutic  agents.  Our  findings 
highlights  the  potential  importance  of  prescreening  patients  for  their  variant  protein  in  order  to 
prognosticate  whether  a  particular  small  molecule  inhibitor  will  be  therapeutically  efficacious. 

Impact  on  society  bevond  science  and  technology; 

Nothing  to  report 


5.  CHANGES/PROBLEMS 
Changes  in  approach: 

None 
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Actual  or  anticipated  problems  or  delays: 

None 

Changes  that  had  significant  impact  on  expenditures: 

None 

Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards, 
and/or  select  agents: 

None 


6.  PRODUCTS 

Publications,  conference  papers,  and  presentations: 

i.  Manuscript  in  review.  Wang  B.-D.,  Ceniccola  K.,  Hwang  S.,  Andrawis  R.,  Horvath  A., 
Freeman  J.A.,  Knapp  S.,  Ching  T.,  Garmire  L.,  Patel  V.,  Garcia-Bianco  M.A.,  Patierno 
S.R.  and  Lee  N.H.  (2016)  Aberrant  alternative  splicing  in  African  American  prostate 
cancer:  novel  driver  of  tumor  aggressiveness  and  drug  resistance.  Nature  Commun. 
(positively  reviewed  by  reviewers;  have  completed  requested  experiments  and 
resubmitted). 

ii.  Presentation  “Alternative  splicing  of  FGFR3  as  a  mechanism  for  prostate  cancer  health 
disparities”  at  the  9th  AACR  Conference  on  The  Science  of  Cancer  Health  Disparities  in 
Racial/Ethnic  Minorities  and  the  Medically  Underserved,  September  25  -  28,  2016. 

hi.  Manuscript  in  review.  Freedman  JA,  Wang  Y,  Xuechan  L,  Liu  H,  Moorman  P,  George 
DJ,  Lee  NH,  Hyslop  T,  Wei  Q  and  Patierno  S  (2016)  Single  nucleotide  polymorphisms  of 
sternness  pathway  genes  predicted  to  regulate  RNA  splicing,  microRNA  and  transcription 
are  associated  with  prostate  cancer  survival.  Clinical  Cancer  Res.  (submitted) 

iv.  Manuscript  in  review.  Wang  Y,  Lreedman  JA,  Liu  H,  Moorman  P,  Hyslop  T,  George  DJ, 
Lee  NH,  Patierno  S  and  Wei  Q  (2016)  Variants  of  sternness-related  genes  predicted  to 
regulate  RNA  splicing  associated  with  racial  disparities  in  susceptibility  to  prostate 
cancer.  Carcinogenesis  (submitted) 

Websitefsi  or  other  Internet  sitefsi: 

None 

Technologies  or  techniques; 

None 

Inventions,  patent  applications,  and/or  licenses: 

i.  Our  provisional  patent  application  has  been  issued.  Application  number:  61/948,218. 

Date  patent  issued:  09/27/16.  Application  title:  Companion  Diagnostics  for  Cancer  and 
Screening  Methods  to  Identify  Companion  Diagnostics  for  Cancer  Based  on  Splicing 
Variants 


Other  Products: 
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Norman  H  Lee,  PhD 
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design  and  statistical  analysis. 
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Bi-Dar  Wang,  PhD 
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xenograft  assays  and  cell-free  protein  purification  and 
analysis. 

Funding  support: 

Partially  supported  by  GWU  bridge  support 

Name: 

Jacqueline  Olender 

Project  Role: 

PhD  graduate  student 

Nearest  person  month  worked: 

12 

Contribution  to  project: 

Contributed  to  the  cloning  of  variant  cDNAs,  in  vitro  and  in 
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i.  Wang  B.-D.,  Ceniccola  K.,  Hwang  S.,  Andrawis  R.,  Horvath  A.,  Freeman  J.A.,  Knapp  S., 
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Abstract 


Despite  the  current  clinical  challenge  of  reducing  prostate  cancer  (PCa)  disparities, 
molecular  mechanisms  contributing  to  race-related  tumor  aggressiveness  have  not 
been  fully  characterized  and  the  RNA  splicing  landscape  of  PCa  across  racial 
populations  has  not  been  explored.  Here,  we  identify  novel  genome -wide,  race-related 
RNA  splicing  events  as  critical  mechanisms  driving  PCa  aggressiveness  and 
therapeutic  resistance  in  African  American  (AA)  men.  AA-enriched  splice  variants  of 
PIK3CD,  FGFR3,  TSC2  and  RASGRP2  contribute  to  greater  oncogenic  potential 
compared  to  corresponding  European  American  (EA)-expressing  variants.  Ectopic 
overexpression  of  the  newly  cloned  AA-enriched  variant,  PIK3CD-S,  in  EA  PCa  cell 
lines  enhances  AKT/mTOR  signaling  and  increases  proliferative  and  invasive 
capacity  in  vitro  and  confers  resistance  to  selective  PI3K8  inhibitor,  CAL-101,  in 
mouse  xenograft  models.  High  PIK3CD-S  expression  in  PCa  specimens  associates 
with  poor  survival.  These  results  highlight  the  potential  of  RNA  splice  variants  to 
serve  as  novel  biomarkers  and  molecular  targets  for  developmental  therapeutics  in 
aggressive  PCa. 


Introduction 


Prostate  cancer  (PCa)  is  the  most  commonly  diagnosed  cancer  and  the  second  leading 
cause  of  cancer  death  among  American  men  Striking  population  disparities  in  PCa 
risk  and  clinical  outcome  have  been  observed  across  racial  and  ethnic  groups. 
Notably,  African  American  (AA)  men  exhibit  1.6-fold  higher  incidence  and  2.4-fold 
higher  mortality  rates  of  PCa  compared  to  European  American  (EA)  men 
Socioeconomic  factors  remain  a  major  component  accounting  for  the  PCa  disparities 
between  AA  and  EA  populations.  However,  higher  mortality  and  recurrence  rates  are 
still  observed  in  AA  PCa  even  after  adjustment  of  the  socioeconomic  factors  between 
AAs  and  EAs  suggesting  that  intrinsic  differences  in  the  biological  characteristics 
and  genomic  composition  of  AA  and  EA  PCa  may  also  play  a  critical  contributing 
role  in  PCa  disparities  among  racial  groups. 

Alternative  splicing  is  a  post-transcriptional  process  that  occurs  during  the 
maturation  of  mRNA  transcripts,  allowing  the  synthesis  of  alternative  mRNA 
transcripts  that  encode  structurally  and  perhaps  functionally  disparate  protein 
isoforms.  Next-generation  DNA  and  RNA  sequencing  (NSG)  data  suggests  that  more 
than  90%  of  human  genes  undergo  alternative  splicing  and  the  resulting  complexity 
in  the  mRNA  transcriptome  explains  how  '-'20,000  human  protein-coding  genes  in  the 
genome  can  lead  to  the  diversity  of  >250,000  distinct  proteins  in  the  proteome. 
Accumulating  evidence  indicates  that  alternative  and/or  aberrant  splicing  of  precursor 
(pre)-mRNA  also  plays  an  important  but  largely  underappreciated  role  in  human 
cancers  including  PCa  pathogenesis  Eor  example,  the  B-cell  lymphoma  2-like  1 
{BCL2L1)  pre-mRNA  is  alternatively  spliced  into  two  variants,  Bcl-xS  and  Bcl-xL, 
encoding  protein  isoforms  with  opposite  biological  effects.  Bcl-xS  is  a  pro-apoptotic 
protein,  while  BcE-xE  has  anti-apoptotic  properties  that  confer  resistance  to 
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chemotherapeutic  agents  in  PCa  eell  line  PC-3  The  fibroblast  growth  faetor 
reeeptor  2  {FGFR2)  pre-mRNA  also  undergoes  alternative  splieing,  where  FGFR2- 
Illb  is  predominately  expressed  in  epithelial  cells  and  FGFR2-IIIc  is  primarily 
assoeiated  with  mesenehymal  eells  and  epithelial-to-mesenehymal  transition  of 
prostate  tumor  cells  Another  example  is  the  TMPRSS2-ERG  gene  fusion 
commonly  found  in  PCa  and  associated  with  poor  clinical  outcome  In  a 
eomparative  study  of  TMPRSS2-ERG  variants  ectopically  overexpressed  in  primary 
prostatie  epithelial  cells,  variants  eontaining  exon  8  were  assoeiated  with  increased 
eell  proliferation,  invasion  and  motility  The  androgen  receptor  (AR)  signaling 
pathway  is  critieally  assoeiated  with  growth  and  survival  in  PCa 
Aberrant/alternative  spliee  variant  is  overexpressed  in  the  hormone-refractory 

PCa,  being  correlated  with  poor  survival  and  higher  reeurrence  rates  following 
surgieal  treatment  Finally,  extensive  unsplieed  pre-mRNAs  have  been  observed  in 
metastatic  castration-resistant  PCa  . 

Despite  the  signifieance  of  alternative  and/or  aberrant  splicing  in  PCa 
development  and  progression  irrespeetive  of  race,  the  occurrence  of  race-specific/- 
enriched  (AA  vs.  EA)  PCa  splicing  events  and  a  causal  relationship  between  these 
events  and  observed  PCa  health  disparities  among  racial  groups  remains  unexplored. 
For  example,  it  is  unclear  if  FGFR2-IIIc,  TMPRSS2-ERG  variants  eontaining  a  72- 
base  exon,  AR-V7  and/or  other  as  yet  undiscovered  variants  assoeiated  with  more 
aggressive  PCa  might  be  predominantly  or  seleetively  expressed  in  AA  PCa,  thus 
possibly  contributing  to  the  PCa  disparities.  In  addition,  it  is  unknown  if  differenees 
in  mRNA  splicing  along  racial/population  lines  occur  in  only  a  limited  number  of 
genes  or  more  globally  aeross  the  transeriptome.  If  the  latter,  it  will  be  important  to 
aseertain  whether  these  genome-wide,  differential  splicing  events  are  over- 


represented  within  speeifie  gene  ontologies  (i.e.  proto-oneogenes,  tumor  suppressor 
genes  and/or  signaling  pathways  playing  roles  in  partieular  biologieal  proeesses). 
Lastly,  if  present,  assessment  of  the  funetional  eonsequenees  of  any  raoe-speeifie  (or  - 
enriehed)  splieing  events  will  provide  eritieal  further  insight  into  the 
genetie/moleeular  meehanisms  underlying  PCa  disparities  among  raeial  groups.  To 
this  end,  we  have  applied  a  funetional  genomies  approaeh  to  address  these  questions 
in  biopsy  and  surgieal  samples  of  AA  and  EA  prostate  eaneer.  We  demonstrate  global 
differential  splieing  events  involving  >2,500  genes  in  prostate  tumors  of  AA  versus 
EA.  The  vast  majority  of  genes  exhibiting  raee-speeifie/-enriehed  differential  splieing 
are  assoeiated  with  eaneer  and  eomputational  analysis  identified  >25  eaneer-related 
pathways  with  an  over-representation  of  differential  splieing  events.  Moreover,  we 
have  evaluated  the  biologieal  signifieanee  of  these  splieing  events  in  a  subset  of 
proto-oneogenes  and  identified  and  eloned  a  novel  AA-enriehed  PIK3CD  short 
variant  {PIK3CD-S)  that  imparts  inereased  proliferative  and  invasive  properties  to 
PCa  eell  lines  eompared  to  an  EA-expressing  long  variant  {PIK3CD-L).  Using 
xenograft  and  metastasis  mouse  models,  we  found  that  PIK3CD-S  eneodes  an  isoform 
that  is  resistant  to  small  moleeule  inhibitors  targeting  PIK3CD,  while  the  variant 
eneoded  by  PIK3CD-L  is  sensitive  to  sueh  inhibitors.  These  results  highlight  an 
opportunity  to  utilize  population  level  differenees  in  tumor  biology  to:  i)  diseover 
novel  spliee  variants  that  will  likely  serve  as  novel  preeision  biomarkers  and/or 
moleeular  targets  for  developmental  therapeuties  against  aggressive  PCa  (partieularly 
in  the  AA  population),  ii)  identify  previously  hidden  spliee  variants  eneoding 
oneogenie  signaling  proteins  resistant  to  small  moleeule  inhibitors,  and  iii)  assimilate 
spliee  variant  information  for  prognostieation  of  eaneer  aggressiveness  and/or 
therapeutie  responsiveness. 


Results 


Exon  array  analysis  reveals  genome-wide  differential  splicing  events  in  AA 
versus  EA  PCa.  A  total  of  35  PCa  (20  AA,  15  EA)  and  35  patient-matched  normal 
prostate  (NP)  specimens  (20  AA,  15  EA)  derived  from  chemo-/hormone-/radiation- 
naiVe  patients  were  interrogated  using  the  Affymetrix  Human  Exon  1 .0  ST  GeneChip 
to  assess  differential  splicing  (DS)  events.  Gleason  scores  of  PCa  specimens  (range  6- 
8)  and  patient  ages  (range  49-81  years)  were  not  significantly  different  between  the 
AA  and  EA  cohorts  (P>0.05  Eisher's  exact  test).  In  AA  PCa  vs.  EA  PCa  and  AA  NP 
vs.  EANP,  the  significant  differentially  expressed  exons  (Figs,  la  and  lb)  could  be 
modeled  using  the  Alternative  Splice  ANOVA  approach  into  2,520  and  2,849  DS 
events,  respectively  (Supplementary  Table  SI,  Supplementary  Fig,  SI).  As 
depicted  in  the  Venn  diagram  (Fig.  Ic),  1,876  genes  (i.e.  2,520  -  644)  exhibited  DS 
events  unique  to  AA  PCa  vs.  EA  PCa,  2,205  differentially  spliced  genes  (i.e.  2,849  - 
644)  were  unique  to  AA  NP  vs.  EA  NP  and  644  DS  events  were  in  common  (i.e.  DS 
events  pre-existing  in  AA  NP  vs.  EA  NP  and  preserved  in  AA  PCa  vs.  EA  PCa). 
Examples  of  genes  with  pre-existing  DS  events  included  PIK3CD,  ITGA4  and  MET, 
while  RASGRP2,  NFl  and  BAKl  are  examples  of  differentially  spliced  genes 
occurring  only  in  AA  PCa  vs.  EA  PCa.  In  EA  PCa  vs.  EA  NP  and  AA  PCa  vs.  AA 
NP,  the  significant  differentially  expressed  exons  (Supplementary  Fig,  S2)  could  be 
modeled  into  1,297  and  1,733  DS  events,  respectively  (Fig.  Ic).  Presumably,  a  subset 
of  1,575  genes  (i.e.  1,733  -  158)  with  DS  events  unique  to  AA  PCa  may  contribute  to 
PCa  disparities.  Examples  in  this  category  included  FGFR3  and  TSC2 
(Supplementary  Table  SI).  On  the  other  hand,  a  subset  of  158  genes  with  DS  events 
in  common  to  both  AA  and  EA  PCa  may  contribute  to  PCa  progression  regardless  of 
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race  (Fig,  Ic).  Consistently,  such  genes  ineluded  TMPRSS2  andAR  (Supplementary 


Table  SIC  and  D).  Analysis  of  the  exon  array  data  employing  both  gene -wise 
and  Alternative  Spliee  ANOVA  modeling  approaehes  identified  898  genes  (i.e. 
1,188  -  290)  that  were  differentially  expressed  but  not  exhibiting  DS  in  AA  PCa  vs. 
EA  PCa,  and  2,230  (i.e.  2,520  -  290)  genes  undergoing  DS  but  not  differential 
expression  (e.g.  level  of  variant  ‘A’  for  gene  ‘X’  in  AA  PCa  equivalent  to  variant  ‘B’ 
for  gene  ‘X’  in  EA  PCa;  Fig,  Ic). 

Prevalence  of  differential  splicing  events  in  cancer-associated  pathways.  We 

eategorized  genes  undergoing  DS  in  AA  PCa  vs.  EA  PCa  based  on  moleeular 
funetion,  gene  ontology  (GO)  and  disease  assoeiation.  Relevant  cancer-related 
ontologies  ineluded  eell  growth  and  proliferation,  eell  death  and  survival,  eellular 
movement,  eell  adhesion  and  DNA  damage/repair  (E-value  ranged  from  6.54x10'^^- 
1.88x10'^;  Supplementary  Table  S2).  Notably,  a  large  fraetion  (1,816  out  of  2,520, 
71.8%)  of  the  differentially  splieed  genes  were  diseovered  to  be  over-represented 
aeross  multiple  eaneers,  ineluding  gastrointestinal,  eoloreetal,  renal,  breast,  brain, 
reetal,  brain,  lung,  stomaeh,  prostate  and  hematologic  cancers  (E-values  ranged  from 
1. 43  xlO'*^- 1.96x10'^;  Fig,  Id  and  Supplementary  Table  S2).  There  was  an 
unexpected  skewing  in  the  distribution  of  in- frame  versus  out-of-frame  exon  skipping 
events  in  eaneer-related  genes,  where  in-frame  events  were  signifieantly  favored  in 
AA  over  EA  PCa  speeimens  (E  <  0.05  Eisher’s  exact  test;  Supplementary  Table  S3). 
This  finding  was  in  line  with  an  overall  signifieant  preferenee  for  in-frame  events 
aeross  all  genes  (eaneer-related  and  non-eancer-related)  in  AA  PCa  specimens 
(Supplementary  Table  S3).  In  the  ease  of  non-eaneer-related  genes  only,  there  was 
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no  significant  skewing  of  in-frame  distribution  events  between  AA  versus  EA  PCa 

(Supplementary  Table  S3). 

We  also  examined  the  distribution  of  DS  events  across  cell  signaling 
pathways.  There  was  a  striking  signifieant  over-representation  of  DS  events  in 
multiple  oncogenic  signaling  pathways,  including  EGF,  VEGE,  PTEN,  PI3K/AKT, 
ERK/MAPK,  cell  cycle  andNExE  signaling  (Eisher’s  Exact  Test,  P-values  of 
0.00126,  0.00891,  0.00426,  0.01288,  0.01622,  0.01259  and  0.02089,  respectively; 
Supplementary  Fig,  S3  and  Supplementary  Table  S4).  Interestingly,  many  of  these 
same  pathways  are  known  to  be  mutated  based  on  earlier  eancer  genome  sequeneing 
studies  A  composite  oneogenie  signaling  pathway  comprised  of  DS  events  found 
in  AA  PCa  vs.  EA  PCa  is  depieted  in  Fig  2.  Taken  together,  our  data  provide  strong 
evidence  that  DS  events  may  play  a  eritieal  role  in  PCa  disparities. 

Validation  of  alternative  splice  variants  in  AA  vs  EA  PCa,  We  proceeded  to 
validate  a  subset  of  both  proto-oneogenes  and  tumor  suppressor  genes  with  DS  events 
in  our  eomposite  cancer  signaling  pathway,  including  PIK3CD,  FGFR3,  TSC2, 
ITGA4,  MET,  NFl,  BAKl,  ATM  and  RASGRP2  (Fig,  2).  RT-PCR  was  performed  on 
RNA  samples  obtained  from  AA  and  EA  PCa  speeimens  originally  interrogated  by 
the  exon  arrays.  Primer  pairs  or  trios  were  designed  for  RT-PCR  to  amplify 
simultaneously  multiple  variants  of  eaeh  gene  (Fig,  3a).  As  shown  in  Fig,  3b,  AA 
PCa  speeimens  eontained  both  PIK3CD  long  {PIK3CD-L,  including  exon  20)  and 
short  {PIK3CD-S,  missing  exon  20)  variants,  whereas  EA  PCa  samples  predominately 
expressed  the  PIK3CD-L  variant.  Hence,  the  RT-PCR  results  were  in  agreement  with 
the  exon  array  data,  indicating  the  presence  of  an  AA-enriched  PIK3CD-S  variant. 
Analogous  findings  were  obtained  where  either  a  short  or  long  variant  of  each  gene 
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was  confirmed  by  RT-PCR  to  be  enriehed  or  uniquely  expressed  in  AA  (TSC2-S, 
ITGA4-L,  MET-L,  BAKl-L)  or  EA  PCa  speeimens  (FGFR3-L,  ITGA4-S,  MET-S, 
NFl-S,  BAKl-S)  (Fig,  3b;  see  Supplementary  Fig,  S4  for  quantitative  RT-PCR 
results  from  n  =  22-25  AA  and  n  =  21-24  EA  PCa  speeimens).  Exon  array  data  also 
revealed  two  alternative  RASGRP2  transeripts  with  apparent  mutually  exelusive  exon 
skipping  events  (Fig,  3a  and  Supplementary  Fig,  Sid).  RT-PCR  validation  likewise 
eonfirmed  that  a  RASGRP2-b  variant  (exon  1 1  exeluded)  was  exelusively  expressed 
in  AA  PCa,  while  a  RASGRP2-a  variant  (exon  12  exeluded)  was  enriehed  in  EA  PCa 
(Fig,  3b  and  Supplementary  Fig,  S4).  We  were  unable  to  validate  differential 
splieing  oiATM,  whereas  2  additional  genes  {GSK3A  and  EPHAl)  identified  by  exon 
arrays  as  not  undergoing  differential  splieing  were  eonfirmed  by  RT-PCR.  In 
summary,  there  was  strong  agreement  (10/11  or  91%)  between  exon  array  and  RT- 
PCR  results,  thus  providing  an  internal  quality  metrie  to  our  global  analysis  of 
differential  splieing  events  between  AA  and  EA  PCa  (see  Supplemental  Materials  & 
Methods  for  additional  metries). 

The  raee-dependent  expression  oiPIK3CD  variants  was  partieularly 
interesting  owing  to  reeent  findings  implieating  PI3K5  (pi  105)  kinase  aetivity  in 
hematologie  malignaneies  as  well  as  other  eaneer  types  In  a  separate  eohort  of 
PCa  speeimens  obtained  from  32  AAs  (age  range  52-76  years,  Gleason  seore  range  6- 
8)  and  30  EAs  (age  range  50-82  years,  Gleason  seore  range  6-8;  not  signifieantly 
different  from  AA,  Eisher’s  exaet  test,  P>0.05),  quantitative  RT-PCR  validation  was 
performed  reaffirming  signifieantly  higher  levels  of  PIK3CD-S  relative  to  PIK3CD-L 
in  AA  eompared  to  EA  PCa  speeimens  (Fig,  3c).  Given  the  robustness  and  potential 
signifieanee  of  these  findings,  subsequent  in  vitro  and  in  vivo  studies  eentered  on  the 
PIK3CD  variants,  as  deseribed  below. 


Molecular  cloning  of  PIK3CD  splice  variants.  The  AA-enriched  PIK3CD-S  variant 
has  never  before  been  deseribed  in  the  literature  nor  the  UCSC  (genome.uesc.edu/)  or 
Ensembl  Genome  Browser  (www.ensembl.org/).  Consequently,  we  cloned  the  full- 
length  version  of  the  PIK3CD-S  variant  from  the  AA  PCa  cell  line  MDA  PCa  2b  and 
PIK3CD-L  variant  from  MDA  PCa  2b  as  well  as  EA  PCa  cell  lines  VCaP  and  ENCaP 
using  standard  molecular  approaches  (e.g.  5’-  and  3 ’-RACE;  see  ^^).  We  likewise 
cloned  matching  PIK3CD-S  and  PIK3CD-L  variants  from  PCa  patient  specimens. 
Supplementary  Fig,  S5  schematically  depicts  the  full-length  clones  of  the  PIK3CD-L 
variant  (comprising  a  total  of  24  exons)  along  with  three  different  AA  PIK3CD-S 
variants  (variant  excluding  exon  8,  variant  excluding  exon  20  and  variant  excluding 
both  exons  8  and  20)  and  one  AA  large  deletion  variant  of  the  PIK3CD  gene. 
Interestingly,  exclusion  of  exon  8  eliminates  a  30  amino  acid  segment  situated 
between  the  Ras-binding  and  C2  domains,  while  exclusion  of  exon  20  deletes  a  56 
amino  acid  segment  located  in  the  catalytic  domain  of  PI3K8.  In  subsequent 
functional  studies  involving  ectopic  overexpression  of  the  short  variant  (see  below), 
we  concentrated  our  efforts  on  the  variant  missing  exon  20  given  the  possibility  that 
kinase  activity  may  be  affected. 

Enrichment  of  PIK3CD-S  isoform  is  associated  with  enhanced  invasive  and 
proliferative  phenotypes  in  AA  PCa  cells.  We  hypothesized  that  the  splice  variants 
specific  or  enriched  in  AA  PCa  may  contribute  to  a  more  aggressive  oncogenic 
phenotype.  To  test  this,  we  designed  both  exon-specific  and  exon  junction-specific 
siRNAs  to  target  either  the  PIK3CD-L  or  PIK3CD-S  variant,  respectively,  in  EA  and 
AA  PCa  cell  lines  and  examined  the  functional  consequences  of  these  knockdowns  on 


both  cell  proliferation  and  invasion.  A  similar  strategy  was  applied  to  investigate  the 
biological  significance  of  the  variants  of  FGFR3,  TSC2  and  RASGRP2.  VCaP  and 
MDA  PCa  2b  cell  lines  were  used  as  population-specific  PCa  models,  as  these  two 
cell  lines  represent  bone  metastases  derived  from  castration-resistant  EA  and  AA  PCa 
patients,  respectively  Transfection  of  VCaP  cells  with  exon  20-specifio  siRNA 
(siP2o)  suecessfully  knoeked-down  PIK3CD-L  expression  by  >8-fold  compared  to 
nonsense  (ns)  siRNA  (see  Fig,  4a,  left  panel),  resulting  in  a  signifieant  loss  of 
proliferative  and  invasive  function  in  VCaP  cells  (Fig,  4b,  left).  Conversely,  in  MDA 
PCa  2b  cells,  a  >5-fold  knockdown  of  PIK3CD-L  increased  the  ratio  of  PIK3CD- 
S/PIK3CD-L  expression  by  nearly  2-fold  (Fig,  4a,  right  panel;  see  1.88  S/L  ratio  for 
ns-  versus  3.46  S/L  ratio  for  siP2o-transfected  eells),  and  this  ‘enriehment’  of  the  AA- 
enriehed  PIK3CD-S  variant  subsequently  enhaneed  proliferation  and  invasion  of  the 
AA  cell  line  (Fig,  4a,  right  panel).  Moreover,  MDA  PCa  2b  cells  exhibited 
significantly  higher  basal  invasive  and  proliferative  capaeities  eompared  to  VCaP 
cells  (see  proliferation  and  invasion  of  siNS-transfected  MDA  PCa  2b  vs.  siNS- 
transfeeted  VCaP;  Fig,  4b,  left  and  right  panels).  To  further  evaluate  the  functional 
impact  of  PIK3CD-S  expression  on  eell  proliferation  and  invasion,  the  EA  and  AA 
PCa  eell  lines  were  transfeeted  with  siPj  (siRNA  speeifieally  targeting  the  junction 
sequence  exons  19  and  21).  Transfection  of  siPj  had  no  effect  on  VCaP  proliferation 
and  invasion,  as  expected  since  this  EA  line  does  not  signifieantly  express  PIK3CD-S 
(Figs,  4a  and  4b,  right  panels).  Upon  the  transfection  of  MDA  PCa  2b  cells  with  siPj, 
PIK3CD-S  expression  was  significantly  knocked-down  (Fig,  4a,  right),  resulting  in  a 
loss  of  cell  proliferation  and  invasion  (Fig,  4b,  right).  Taken  together,  these  results 
suggest  that  PIK3CD-S  is  the  more  aggressive  variant,  promoting  PCa  proliferation 
and  invasion  to  a  greater  extent  than  PIK3CD-L. 


Several  additional  exon-speeifie  siRNAs  were  designed  to  test  whether  other 
AA-speeilie/-enriehed  spliee  variants  also  funetionally  eontribute  to  greater  PCa 
aggressiveness.  SiRNAs  targeting  exon  14  {siFGFR3-exl4),  exon  20  {siTSC2-ex20) 
and  exon  1 1  {siRASGRP2-exl  1)  were  used  to  seleetively  suppress  expression  of 
FGFR3-L,  TSC2-L  and  RASGRP2-a  spliee  variants  (predominately  expressed  in  EA), 
respeetively.  Upon  siRNA-mediated  knoekdown  of  FGFR3-L,  TSC2-L  or  RASGRP2- 
a  (Supplementary  Figs,  S6a-c,  top  panels)  in  MDA  PCa  2b  eells,  the  expression 
ratios  FGFR3-S/FGFR-L,  TSC2-S/TSC2-L  aad  RASGRP2-b/RASGRP2-a  inereased 
and  eorrelated  with  augmented  invasive  and/or  proliferative  eapaeity  of  the  AA- 
derived  MDA  PCa  2b  eells  (Supplementary  Figs.  S6a-c,  bottom  panels). 

Colleetively,  our  in-vitro  studies  strongly  suggest  that  the  AA-enriehed  spliee  variants 
PIK3CD-S,  FGFR3-S,  TSC2-S  and  RASGRP2-b  promote  PCa  aggressiveness. 

PIK3CD-S  isoform  promotes  activation  of  AKT/mTOR  signaling.  As  PI3K  plays  a 
central  role  in  the  PI3K/ AKT/mTOR  signaling  pathway,  we  examined  the  ability  of 
the  alternative  PI3K6  isoforms  (encoded  by  the  PIK3CD-L  and  PIK3CD-S  splice 
variants)  to  activate  down-stream  signaling  components  within  this  pathway.  SiRNA 
(siP2o)-mediated  knockdown  of  PIK3K-L  expression  (confirmed  by  qRT-PCR)  led  to 
a  drastic  decrease  in  AKT  phosphorylation  at  Thr308  and  Ser473  while  moderately 
decreasing  (-'2-fold)  phosphorylation  of  mTOR  in  EA  cell  line  VCaP  (Fig,  4c,  top 
panel).  In  contrast,  knockdown  of  PIK3CD-L  in  AA  MDA  PCa  2b  cells  (resulting  in 
an  increased  PIK3CD-S/-L  ratio  confirmed  by  qRT-PCR)  led  to  a  sizable  increase  (2- 
3-fold)  in  phosphorylation  of  AKT,  mTOR  and  S6  (Fig,  4c,  top  panel). 

In  parallel  experiments,  siRNA  (siPj)-mediated  knockdown  of  PIK3K-S 
expression  (confirmed  by  qRT-PCR)  in  MDA  PCa  2b  cells  resulted  in  a  drastic 
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decrease  in  the  phosphorylation  status  of  AKT,  mTOR  and  ribosomal  protein  S6  (Fig, 
4c,  bottom  panel).  As  expected,  treatment  of  VCaP  cells  with  the  siRNA  siPj  had 
negligible  effects  on  AKT,  mTOR  and  ribosomal  protein  S6  phosphorylation,  as  this 
EA  line  does  not  significantly  express  PIK3CD-S.  Taken  together,  the  distinct 
phosphorylation  patterns  of  AKT,  mTOR  and  ribosomal  protein  S6  in  AA  and  EA 
PCa  cell  lines  upon  selective  knockdown  of  either  PIK3CD-L  or  PIK3CD-S,  again 
suggested  that  PIK3CD-S  is  the  more  aggressive  variant,  promoting  oncogenic 
signaling. 

PI3K8-L  and  PI3K8-S  isoforms  exhibit  differential  response  to  small  molecule 
inhibitors  in  vitro  and  in  mouse  xenograft  and  metastasis  models.  We  tested 
whether  pharmacological  inhibition  of  PI3K8  isoforms  represented  a  potential 
strategy  for  ameliorating  PCa  aggressiveness.  CAL-101,  a  small  molecule  inhibitor 
(SMI)  specific  for  PI3K8  was  employed  to  assess  its  inhibitory  effects  on 
oncogenic  signaling  and  proliferation  in  EA  PCa  cell  lines  VCaP  and  PC-3  that  stably 
overexpressed  the  His-tagged  PI3K5-S  (excluding  exon  20)  or  PI3K5-E  isoform 
(including  exon  20).  Equivalent  levels  of  PI3K5  isoform  expression  in  each  cell  line 
was  confirmed  by  western  blot  with  a  His  tag  antibody  (Fig,  5a).  In  both  EA  PCa  cell 
lines,  ectopic  overexpression  of  PI3K5-S  was  associated  with  a  2-3  fold  greater 
phosphorylation  of  AKT  and  ribosomal  protein  S6  compared  to  ectopic 
overexpression  of  PI3K6-E  (Fig,  5a,  absence  of  CAL-101  treatment).  CAL-101 
(50mg/kg)  induced  a  significant  reduction  in  basal  AKT  and  ribosomal  protein  S6 
phosphorylation  (Fig.  5a)  and  a  dose-dependent  inhibition  of  proliferation  in  both  EA 
cell  lines  overexpressing  the  PI3K5-L  variant  (Fig.  5b).  In  contrast,  CAL-IOI 
(50mg/kg)  had  negligible  effects  on  inhibiting  basal  AKT  signaling  in  EA  PCa  cell 


lines  overexpressing  PI3K5-S,  as  phosphorylation  states  of  AKT  and  ribosomal 
protein  S6  were  eomparable  to  vehiele  treated  eells  (Fig.  5a).  In  agreement,  BrdU 
labeling  assays  demonstrated  that  proliferation  of  VCaP  and  PC-3  eells  eetopieally 
overexpressing  PI3K6-S  was  greater  than  eells  overexpressing  PI3K5-L 
(Supplementary  Fig.  S7).  Moreover,  PI3K5-S-overexpressing  VCaP  and  PC-3  eells 
were  not  effeetively  inhibited  by  CAL-101  treatment;  only  at  extreme  doses  of  CAL- 
101  (>30pM)  was  proliferative  aetivity  in  PC-3  eells  signifieantly  impaired  (Fig.  5b). 
In  eontrast,  the  AKT  inhibitor  MK-2206  dose-dependently  deereased  proliferation 
in  both  PI3K6-S-  and  PI3K5-L-overexpressing  VCaP  and  PC-3  eells  (Fig.  5c).  These 
results  suggest  that  PI3K5-S-stimulated  proliferation  is  resistant  to  CAL-101 
inhibition  in  sharp  eontrast  to  PI3K5-L;  while  inhibition  of  AKT,  whieh  is 
downstream  of  PI3K5-S,  effeetively  bloeked  proliferation. 

To  examine  the  effeets  of  PIK3CD  spliee  variants  on  tumor  growth  in  vivo,  we 
subeutaneously  injeeted  2x10®  PC-3  eells  stably  overexpressing  equivalent  amounts 
of  the  PI3K5-L  or  PI3K6-S  (missing  exon  20)  isoform  into  the  left  hind  flank  of 
NOD-SCID  miee.  Miee  harboring  PI3K5-L-expressing  or  PI3K5-S-expressing  PC-3 
eell  xenografts  were  administrated  either  vehiele  (PBS)  or  CAL-101  (50  mg/kg)  by 
daily  i.p.  injeetion.  CAL-101  treatment  for  30  days  signifieantly  redueed  the  growth 
of  PI3K5-L-expressing  xenografts  eompared  to  the  vehiele  treatment  (Figs.  6a  and 
6b).  In  eontrast,  miee  with  xenografts  of  PI3K5-S-expressing  eells  had  negligible 
suppression  of  their  xenograft  growth  following  CAL-101  treatment  eompared  to 
vehiele-treated  animals  (Figs.  6a  and  6b). 

We  further  examined  the  inhibitory  effeets  of  CAL-101  on  PI3K6  isoforms  in 
an  in  vivo  tumor  metastasis  model.  1x10^  PIK3CD-L-  or  T’/KiCD-S'-overexpressing 
PC-3  eells  were  injeeted  into  the  tail-vein  of  NOD-SCID  miee,  and  animals  were 


subsequently  administrated  with  vehiele  or  CAL-101  (50mg/kg)  via  i.p.  injeetion  (3 
times  per  week).  After  8  weeks,  vehiele -treated  miee  earrying  PI3K5-L- 
overexpressing  eells  developed  prominent  tumor  metastases  in  the  lungs  (Fig,  6c  and 
6d),  while  the  CAL-101  treatment  group  exhibited  a  >50%  reduetion  {P  <  0.05)  of 
metastases  (Figs,  6c  and  6d).  By  comparison,  CAL-101  treatment  failed  to 
significantly  inhibit  tumor  metastases  in  mice  harboring  PI3K5-S-expressing  cells 
(Figs,  6c  and  6d).  Noteworthy,  the  size  of  lung  metastases  (average  area  of  nodules) 
in  mice  harboring  PI3K5-S-overexpressing  cells  was  slightly  greater  (-'15%) 
compared  to  animals  with  PI3K5-L-overexpressing  cells,  albeit  not  statistically 
significant  (P>0.05).  Taken  together,  the  in  vitro  and  in  vivo  functional  studies 
suggest  that  SMIs  such  as  CAL-101  (i.e.  competitive  ATP  binding  inhibitors 
may  be  ineffective  against  the  PI3K5-S  isoform  in  AA  PCa. 

Cell-free  PI3K6  isoform  kinase  assay.  The  result  of  the  exclusion  of  exon  20  (168 
bp)  in  the  PIK3CD-S  variant  is  an  in-frame  deletion  of  56  amino  acids  (residues  810- 
865)  in  the  catalytic  domain  of  the  PI3K6-S  isoform  (Fig,  7a).  To  gain  further  insight 
into  the  functional  differences  between  PI3K5  isoforms,  the  interaction  of  PI3K5-L 
and  -S  with  regulatory  subunit  p85a  was  investigated.  Whole  cell  lysates  from  PC-3 
cells  overexpressing  either  His-tagged  PI3K5-S  or  PI3K5-L  were  subjected  to  western 
analysis,  demonstrating  that  each  cell  line  expressed  equivalent  levels  of  their 
respective  PI3K6  isoform  as  well  as  equal  p85a  expression  (Fig,  7b;  left  panel). 
Interestingly,  co-IP  of  the  PI3K6/p85a  complex  from  whole  cell  lysates  using  an  anti- 
His  antibody  demonstrated  that  p85a  bound  with  3-4-fold  greater  proficiency  to 
PI3K5-L  compared  to  p85a  binding  to  PI3K5-S  (Fig,  7b;  right  panel,  column  E). 


Binding  proficiency  was  inversely  correlated  with  PI3K6  isoform  kinase  activity  (Fig. 
7c;  right  panel). 

Next,  His-tagged  PI3K5  isoforms  were  purified  from  the  lysates  of  PC-3  cells 
overexpressing  either  PI3K5-S  or  PI3K5-L  using  Ni-NTA  resin  columns.  As  shown  in 
Fig,  7d  (left  and  middle  panels),  purification  of  the  PI3K5-S  and  -L  isoforms  was 
verified  by  western  blotting  using  anti-His  or  anti-PI3K6  antibodies.  In  addition,  the 
Ni-NTA  resin  column  purification  approach  resulted  in  the  isolation  of  PI3K6 
isoforms  that  were  no  longer  bound  to  p85a  (Fig.  7d;  far  right  panel,  column  E).  The 
purified  PI3K6  isoforms  (minus  p85a)  were  incubated  with  vehicle,  non-selective 
PI3K  inhibitor  wortmannin  (100  nM)  or  PI3K6-specific  inhibitor  CAL-101  (100 
nM),  and  subjected  to  a  PI3K  activity  assay.  In  the  absence  of  bound  p85a,  kinase 
activity  of  PI3K5-L  was  equivalent  to  PI3K5-S  (Fig,  7e,  compare  vehicle  treatments). 
In  agreement,  siRNA-mediated  knockdown  of  p85a  in  wild-type  EA  PCa  cell  lines 
VCaP  and  PC-3  was  associated  with  an  increase  in  invasive  activity  (Supplementary 
Fig,  S8).  Remarkably,  wortmannin  and  CAL-101  significantly  inhibited  the  activity 
of  the  PI3K6-L  isoform,  but  not  the  PI3K5-S  isoform  (Fig,  7e).  These  results 
demonstrate  that  PI3K5-S  maintains  kinase  activity  even  in  the  presence  of  SMIs, 
supporting  the  in-vitro  and  in-vivo  results  shown  in  Figs,  5  and  6. 

Discussion 

Emerging  evidence  indicates  that,  in  addition  to  socioeconomic  factors 
biologic/genetic  factors  likely  contribute  to  PCa  health  disparities  among  racial 
groups  Immunohistochemical  evaluation  has  revealed  that  epidermal  growth 
factor  receptor  (EEGR)  overexpression  in  PCa  is  more  common  in  AA  (-'50%  of 
patients)  compared  to  EA  patients  (<20%)  Moreover,  a  recent  study  in  our 


laboratory  demonstrated  that  EGFR  signaling  is  exaeerbated  by  the  presenee  of  AA- 
speeifie  miRNA-mRNA  interaetions,  whieh  appears  to  promote  PCa  aggressiveness 
in  AA  patients  .  In  addition  to  overexpression  of  EGFR,  androgen  reeeptor  (AJi) 
mRNA  levels  are  81%  higher  in  AA  PCa  eompared  to  EA  PCa  .  Interestingly,  a 
polymorphie  CAG  repeat  in  exon  1  of  the  AR  gene  is  frequently  observed  to  be 
shorter  in  AA  PCa  eompared  to  eorresponding  EA  speeimens  The  shorter  repeats 
are  assoeiated  with  higher  levels  of  AR  gene  transeription  and  poorer  prognostie 
indiees,  sueh  as  higher  eaneer  grade  and  stage,  along  with  inereased  ineidents  of 
metastasis  and  mortality  .  Eastly,  several  gene  profiling  studies  have  identified 
eytokine,  immune  response,  inflammatory,  ehemotaxis  and  pro-metastasis  signaling 
pathways  that  are  signifieantly  enriehed  with  differentially  expressed  genes  in  AA 
PCa  eompared  to  EA  PCa 

In  addition  to  DNA  mutation  and  gene  expression  ehanges,  alternative  and/or 
aberrant  RNA  splieing  has  been  demonstrated  to  augment  oneogenie  aetivity  (e.g.  eell 
proliferation,  invasion,  motility  and  therapeutie  resistanee)  Alternative  splieing 
of  the  gene  eneoding  the  Bel-x  apoptotie  regulator  results  in  anti-apoptotie  Bel-XL  and 
pro-apoptotie  Bel-xs  isoforms  PCa  xenografts  expressing  Bel-XL  have  been  shown 
to  exhibit  eastration-resistant  tumor  growth,  and  high  levels  of  Bel-XL  are  assoeiated 
with  metastatie,  eastration-resistant  human  tumor  samples  Manipulation  of  splieing 
to  deerease  Bo1-xl  and  inerease  Bel-xs  levels  has  been  shown  to  reduee  tumor  load  in 
a  murine  melanoma  tumor  model  Aberrant  splieing  of  the  androgen  reeeptor  pre- 
mRNA  generates  variants  (e.g.  AR-V7)  laeking  the  C-terminal  ligand  binding  domain 
AR-V7  has  been  assoeiated  with  eastrate  resistant  PCa,  poorer  elinieal  outeomes 
and  resistanee  to  androgen  ablation/androgen  reeeptor  inhibition  therapies 


The  phenomenon  of  DS,  mueh  less  global  DS  events,  has  not  been  adequately 
explored  as  a  possible  meehanism  underlying  PCa  health  disparities  among  raeial 
groups.  Potential  involvement  of  the  constitutively  aetive  AR-V7  spliee  variant  in  PCa 
health  disparities  has  been  suggested  in  a  reeent  study.  Seleetive  down-regulation  of 
miR-212  observed  in  AA  PCa  eompared  to  EA  PCa  is  correlated  with  up-regulation 
of  splicing  factor  hnRNP-H  1  protein  levels  (miR-212  is  predicted  to  target  hnRNP-Hl 
mRNA),  up-regulation  of  the  AR-V7  variant,  and  resistance  to  anti-androgen  therapy 
in  PCa  cell  lines  In  contrast  to  this  localized  splicing  event,  our  study  reveals  that 
DS  on  a  global  scale  may  be  a  critical  molecular  mechanism  underlying  PCa  health 
disparities.  In  a  comparison  of  AA  PCa  vs.  EA  PCa,  DS  events  were  found  to  be 
highly  prevalent  in  cancer-associated  genes  and  pathways  (Supplementary  Table  SI 
and  S2).  Interestingly,  the  number  of  genes  harboring  DS  events  (2,520  genes)  was  ~3 
times  greater  than  the  number  of  differentially  expressed,  but  not  differentially 
spliced,  genes  (886  genes).  These  findings  have  two  major  implications.  Eirst, 
alternative  and/or  aberrant  splicing  of  pre-mRNAs  may  have  a  greater  role  than 
differential  gene  expression  in  driving  PCa  health  disparities.  Second,  DS  events 
identified  in  our  study  were  statistically  over-represented  in  oncogenic  signaling 
pathways.  In  many  cases,  these  same  pathways  are  known  to  harbor  a  preponderance 
of  gene  mutations  across  different  cancer  types  Thus,  future  studies  are 
warranted  to  investigate  the  oncogenic  consequences  of  interactions  between  gene 
mutations  and  DS  events  within  the  same  pathway.  Taken  together,  DS  adds  another 
layer  of  complexity  to  the  existing  molecular  repertoire  (e.g.  differential  gene 
mutation,  expression,  methylation;  see  review  '^^)  driving  AA  PCa  aggressiveness. 

Studies  on  alternative  splicing  indicate  that  approximately  half  of  such  events 
occurring  in  the  coding  sequence  are  in-frame,  while  the  remaining  events  are  frame- 


shifts  leading  to  truncated  or  extended  C-terminal  proteins  Remarkably,  70%  of 
AA-enriched/-specific  variants  in  our  composite  oncogenic  signaling  pathway  (Fig.  2) 
were  in-frame,  including  PIK3CD-S,  FGFR2-S,  FGFR3-S,  TSC2-S,  RASGRP2-b, 
ATMS  and  GSK3S  (see  Supplementary  Table  S5).  In  stark  contrast,  only  27.3  %  of 
EA-enriched/-specilic  DS  events  in  our  composite  oncogenic  signal  pathway 
exhibited  in- frame  preservation  (Supplementary  Table  S5),  while  the  remaining 
72.7%  of  EA-enriched  DS  variants,  including  ITGA4S,  METS,  NFIS,  RASGRP2-a, 
mTOR-S  and  BAKl-S,  were  frame-shifted.  Why  the  vast  majority  of  DS  events  appear 
to  be  in- frame  in  AA  PCa  vs.  frame-shifted  in  EA  PCa  remains  unresolved 
(Supplementary  Table  S3).  Presumably,  the  preponderance  of  AA  in-frame  events 
detected  in  oncogenic  signaling  pathways  may  be  contributing  to  the  more  aggressive 
nature  of  AA  PCa.  Possible  mechanisms  that  could  drive  differences  in  alternative 
splicing  events  include  differential  expression  of  trans-acting  splicing  factors 
and/or  single  nucleotide  polymorphisms  in  cA-acting  splicing  elements  of 
alternatively  spliced  genes  In  fact,  a  number  of  splicing  factor  mRNAs  appear  to 
be  overexpressed  (e.g.  SRSF2,  SRSF7)  in  AA  PCa  compared  to  EA  PCa 
Regarding  the  in-frame  variants  (i.e.  PIK3CD-L,  FGFR3-L  and  TSC2-L)  detected  in 
EA  PCa  and  tested  experimentally  in  vitro  and/or  in  vivo,  it  is  noteworthy  that  each 
conferred  a  less  aggressive  oncogenic  phenotype  compared  to  the  corresponding  in¬ 
frame  variants  detected  in  AA  PCa  (i.e.  PIK3CD-S,  FGFR3-S  and  TSC2-S). 

Approximately  one-third  of  the  AA-enriched/-specilic  variants  identified  in 
AA  PCa  were  likewise  present  in  patient-matched  NP  specimens,  whereas  the 
remaining  AA-enriched/-specilic  variants  found  in  PCa  were  absent  in  patient- 
matched  NP  specimens  and  thus  appear  to  be  de  novo  events  (i.e.  occurring  as  NP 
evolved  into  PCa).  Accordingly,  the  AA-enriched/-specific  variants  already  present  in 
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NP  specimens  have  the  potential  to  serve  as  inherent  ‘at-risk  alleles’  for  poor  PCa 
prognosis  in  the  AA  population.  By  comparison,  the  de  novo  appearance  of  tumor- 
specific  variants  may  drive  poorer  outcomes.  The  PIK3CD-S  variant  would  be  an 
example  of  a  potential  AA  ‘at-risk  allele’  contributing  to  increased  PCa 
aggressiveness  upon  disease  presentation.  Indeed,  ectopic  overexpression  of  the  AA- 
enriched  PIK3CD-S  variant  in  PCa  cell  lines  was  demonstrated  to  enhance  oncogenic 
potential  (i.e.  increased  invasion,  proliferation  and  AKT/mTOR  signaling)  compared 
to  the  corresponding  EA-enriched  variant  PIK3CD-L.  Moreover,  genetic 
manipulation  of  AA  cell  line  MDA  PCa  2b  PCa  to  favor  expression  of  the  -S  variant 
over  the  -L  variant  likewise  increased  oncogenic  behavior.  Conversely,  genetic 
manipulation  in  the  opposite  direction  had  the  effect  of  decreasing  oncogenic 
behavior  in  the  AA  PCa  cell  line.  Interestingly,  survival  plots  generated  from  TCGA 
RNA-Seq  data  demonstrate  that  a  high  S/L  ratio  is  associated  with  significantly  worse 
survival  for  PCa  and  trending  for  worse  survival  in  both  breast  and  colon  cancer 
(Supplementary  Fig,  S9).  Survival  plots  were  not  stratified  by  race  as  this 
information  is  not  currently  available  in  TCGA.  Given  the  number  patients  analyzed, 
it  seems  highly  probable  that  high  S/L  ratio  values  (besides  being  prevalent  in  AA 
patients)  may  also  be  associated  with  a  subset  of  LA  patients,  suggesting  that 
PIK3CD-S  may  be  useful  in  predicting  survival  in  all  patients  irrespective  of  race. 
Besides  PIK3CD-S,  we  have  identified  an  additional  732  potential  ‘at-risk  alleles’ 

(e.g.  ITGA4-L  and  MET-L)  that  may  be  associated  with  poor  PCa  prognosis  in  the  AA 
population.  Further  experimentation  is  needed  to  investigate  whether  these  variants 
can  serve  as  novel  precision  biomarkers  to  address  PCa  health  disparities  among 
racial  groups.  In  contrast  to  the  ‘at-risk  alleles’,  AA-enriched  variants  FGFR-S  and 
TSC2-S  were  detected  in  AA  PCa,  but  not  in  patient-matched  NP  specimens.  The 
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appearance  of  these  de  novo  variants  during  PCa  formation  may  contribute  to  driving 
the  more  aggressive  PCa  phenotype  observed  in  the  AA  population,  since  in  vitro 
genetic  manipulation  favoring  expression  of  these  -S  variants  over  the  -L  variants 
promoted  oncogenesis  in  AA  cell  line  MDA  PCa  2b.  It  is  noteworthy  that  three  PCa- 
associated  splice  variants  identified  in  previous  studies,  Bcl-xL,  FGFR2-IIIc  and 
TMPRSS2-ERG  variants  containing  exon  8  did  not  exhibit  DS  in  our  AA  PCa 

vs.  EA  PCa  comparison  (see  Supplementary  Table  SI),  suggesting  that  these 
variants  may  contribute  to  PCa  progression  and/or  aggressiveness  in  a  race- 
independent  manner. 

The  identification  of  PIK3CD-S,  a  variant  newly  discovered  and  cloned  in  our 
study,  as  an  ‘at-risk  allele’  for  PCa  aggressiveness  is  germane  given  that  PI3K 
signaling  is  aberrantly  activated  in  a  variety  of  cancers  and  PI3K  inhibitors  have  been 
developed  as  targeted  therapeutics  Class  lA  PBKs  consist  of  three  isoforms, 
including  PI3Ka,  PI3KP  and  PI3K5.  Unlike  PI3Ka  and  PI3KP  that  are  ubiquitously 
expressed,  PI3K6  appears  to  be  preferentially  expressed  in  leukocytes  Previous 
studies  have  revealed  a  crucial  role  of  PI3K5  in  development  and  progression  of 
lymphoid  and  myeloid  malignancies  Interestingly,  accumulating  evidence 
suggests  a  functional  role  of  PI3K5  in  promoting  non-hematologic  tumors  as  well.  For 
example,  overexpression  of  PIK3CD  mRNA  and/or  PBKS  protein  has  been  detected 
in  glioblastoma  ,  neuroblastoma  ,  breast  cancer  and  PCa  ,  and  overexpression 
of  PIK3CD  mRNA  has  been  implicated  in  promoting  cell  growth  and  survival  in 
breast  cancer  and  neuroblastoma  Consistent  with  the  aforementioned  findings, 
our  IHC  experiments  using  a  pan-PBKS  antibody  likewise  revealed  strong  expression 
of  PBKS  protein  in  PCa  specimens  as  well  as  PCa,  breast  cancer  and  colon  cancer 
cell  lines  (Supplementary  Fig.  SIO).  Importantly,  our  study  provides  greater 
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granularity  by  being  the  first  to  demonstrate  the  relationship  between  expression  of  a 
raee-enriehed  PIK3CD  spliee  variant  and  cancer  aggressiveness  as  well  as  resistance 
to  small  molecule  inhibitors  targeting  PI3K6. 

Aberrant  pre-mRNA  splicing  has  recently  been  demonstrated  to  mediate 
therapeutic  resistance  in  multiple  cancer  types.  For  example,  the  constitutively  active 
AR-V7  splice  variant  (lacking  exonic  sequences  encoding  the  ligand  binding  domain) 
confers  resistance  to  enzalutamide  and  abiraterone  acetate  in  castration-resistant  PCa 
patients  In  addition,  melanoma  patients  harboring  BRAF  splice  variants  encoding 
protein  isoforms  that  are  missing  the  RAS-binding  domain  exhibited  resistance  to  the 
RAF  inhibitor  vemurafenib  Noteworthy,  these  studies  did  not  investigate 
whether  variant  expression  and  therapeutic  responsiveness  could  be  stratified  along 
racial  lines.  We  now  provide  evidence  that  the  AA-enriched  variant  PIK3CD-S, 
cloned  in  this  study  for  the  first  time,  imparts  PCa  cell  lines  with  significant  resistance 
to  small  molecule  inhibitors  targeting  PI3K5,  as  demonstrated  in  both  in  vitro  assays 
and  preclinical  mouse  models  of  PCa.  This  short  variant  is  missing  exon  20,  encoding 
a  56  amino  acid  segment  that  is  present  in  PIK3CD-L.  Amino  acids  residing  in  the 
exon  20-encoded  cassette,  appear  to  be  critical  for  the  docking  of  CAL-101  and 
wortmannin.  Indeed,  molecular  modeling  studies  predict  that  Glu826  and  Val828 
(missing  in  PI3K5-S  isoform)  undergo  hydrogen  bonding  with  CAL-101 
Noteworthy,  overall  response  of  indolent  lymphoma  and  chronic  lymphocytic 
leukemia  to  CAL-101  ranges  from  48%  to  81%  Given  our  findings,  it  would  be 
of  interest  to  determine  if  patients  with  primary  resistance  harbor  malignant  cells 
expressing  the  CAL-lOl-resistant  PIK3CD-S  variant,  while  responsive  patients  harbor 
malignant  cells  expressing  the  CAL-lOl-sensitive  PIK3CD-L  variant. 


P85  regulatory  subunits  are  known  binding  partners  of  elass  I  PBKs,  resulting 
in  protein  stabilization  and  inaetivation  of  basal  kinase  aetivity  Somatie 
mutations  in  PIK3R1  (eneoding  p85a)  have  been  identified  that  abrogate  the 
inhibitory  aetion  of  p85a  on  PBKa  in  eaneers  Our  eell-free  assays  demonstrated 
that  p85a  binds  more  effieiently  with  PBK5-L  eompared  to  PBK5-S.  This  interaetion 
appears  to  be  responsible  for  the  lower  kinase  aetivity  exhibited  by  the  long  isoform, 
as  disruption  of  binding  led  to  a  long  isoform  with  inereased  kinase  aetivity 
eomparable  to  the  short  isoform  (see  Fig.  7c  and  e).  The  amino  aeid  Asn334  loeated 
on  the  N-terminal  side  of  PBK6  has  been  postulated  to  serve  as  a  eritieal  eontaet 
point  with  p85a  Our  findings  suggest  that  amino  aeids  810-865,  eneoded  by  exon 
20  and  missing  in  PBK6-S,  may  also  eontain  essential  amino  aeids  required  for 
effieient  eoupling  to  p85a.  Alternatively,  the  presenee  of  amino  aeids  810-865 
permits  PBK6-L  to  adopt  a  eonformation  where  Asn334  (and  other  amino  aeids)  is 
available  to  interaet  with  p85a. 

The  identifieation  and  funetional  validation  of  global  alternative  and  aberrant 
mRNA  splieing  in  eaneer  pathogenesis  remains  ehallenging  and  thus,  largely 
unexplored.  Here,  we  have  undertaken  sueh  an  analysis  in  the  eontext  of  raee-related 
aggressive  PCa  and  have  identified  a  large  number  of  DS  events  in  eaneer-assoeiated 
pathways  in  EA  and  AA  PCa,  with  a  subset  of  these  events  also  being  deteeted  in 
patient-matehed  NP  speeimens.  These  events  will  have  both  biologieal  and  olinieal 
eonsequenees,  ease  in  point  PIK3CD-S.  The  identifieation  of  novel  spliee  variants  as 
biomarkers  and/or  development  of  therapeuties  targeting  protein  isoforms  have  the 
potential  to  reduee  eaneer  health  disparities. 


Methods 


Materials.  EA  PCa  cell  lines  LNCaP  (CRL-1740),  VCaP  (CRL-2876)  and  PC-3 
(CRL-1435),  and  AA  PCa  cell  line  MDA  PCa  2b  (CRL-2422)  were  obtained  from  the 
American  Type  Tissue  Collection  (Manassas,  VA).  Primer  sequences  for  RT-PCR  are 
provided  in  Supplementary  Table  S6.  SiRNAs  were  purchased  from  GE  Dharmaeon 
(Eafayette,  CO)  and  sequences  are  as  follows:  Nonsense,  CCA  AAU  UAU  ACC 
UAC  AUU  GCU;  siP2o,  CCA  ACA  UCC  AAC  UCA  ACA  A;  siPj,  UGA  GGG  AGG 
CCC  UGG  AUC  GA;  siE,  CUC  GAC  UAC  UAC  AAG  AAG  A;  siTSC2-ex20,  CUG 
CGC  UAU  AAA  GUG  CUC  A;  siRASGRP2-exl  1,  CCA  CAU  CUC  ACA  GGA 
AGA  A.  siPIK3Rl  Smart  Pool  (AGU  AAA  GCA  UUG  UGU  CAU  A,  CCA  ACA 
ACG  GUA  UGA  AUA  A,  GAC  GAG  AGA  CCA  AUA  CUU  G,  UAU  UGA  AGC 
UGU  AGG  GAA  A). 

Collection  of  PCa  clinical  specimens.  Prostate  biopsy  samples  were  collected  at  the 
George  Washington  University  Medical  Eaculty  Associates  according  to  an  IRB- 
approved  protocol  (IRB#020867).  High-quality  PCa  and  patient-matched  normal 
prostate  (NP)  biopsy  cores  from  each  of  20  AA  and  15  EA  primary  PCa  patients  were 
collected  and  processed  for  the  exon  array  analysis.  PCa  cores  were  determined  by  a 
pathologist  to  have  Gleason  scores  of  6-7  (17  AA  and  13  EA)  or  8-9  (3  AA  and  2 
EA),  while  NP  cores  were  diagnosed  negative  for  cancer.  There  were  no  significant 
differences  (t-test,  E>0.05)  between  the  two  racial  groups  with  respect  to  age  (average 
age  for  AAs  was  62.3  ±  8.2,  average  age  for  EAs  was  63.3  ±  9.2)  and  Gleason  score 
(range  6-8;  Eisher's  exact  test,  E’>0.05).  No  distant  metastasis  was  detected  in  the 
enrolled  patients. 


Exon  array  and  statistical  analyses.  Total  RNA  was  purified  from  PCa  and  patient- 
matched  NP  biopsy  cores  using  the  RNeasy  micro  kit  (Qiagen,  Valencia,  CA).  High 
quality  RNA  isolation  was  confirmed  by  using  the  Agilent  Bioanalyzer  as  per 
manufacturer’s  protocol  (Agilent  Technologies,  Santa  Clara,  CA).  For  exon  array 
analysis,  1  pg  of  purified  RNA  sample  from  each  biopsy  core  was  interrogated  with 
the  Affymetrix  Human  Exon  1.0  ST  GeneChip  (Santa  Clara,  CA).  Exon  microarray 
data  can  be  assessed  at  GEO  using  accession  number  GSE64331.  The  exon  array  raw 
data  were  subjected  to  quantile  normalization,  GC-content  adjustment,  RMA 
background  correction  and  log2-transformation.  Data  analysis,  data  visualization  and 
statistical  analysis  were  performed  using  Partek  Genomics  Suite  6.6  software  (Partek 
Incorporated,  St.  Eouis,  MO),  as  previously  described  Detection  of  differential 
expression  at  the  gene  level  (gene-wise  analysis)  was  performed  in  Partek  using  a 
One-Step  Tukey’s  Biweight  algorithm  for  detection  of  outlier  probe-sets.  Statistical 
analysis  of  exon  expression  data  was  based  on  ANOVA  with  multiple-correction 
testing  using  10%  Ealse  Discovery  Rate  (EDR)  criterion  DS  events  were  modeled 
using  the  Alternative  Splice  ANOVA  algorithm  implemented  in  Partek  together  with 
selection  of  probe-sets  exhibiting  significant  alternative  splicing  score  determined  at  a 
2%  EDR  Principal  component  analysis  (PCA)  plots  and  two-dimensional  (2D) 
hierarchical  clustering  of  exon  level  data  were  performed  using  Partek.  DS  events 
were  tested  for  statistical  over-representation  in  canonical  signaling  pathways  by 
Eisher’s  exact  test  using  the  Ingenuity  Pathway  Analysis  (IP A)  program  (Ingenuity 
Systems,  Redwood  City,  CA). 

RT-PCR  validation  of  alternative  splice  variants  in  AA  and  EA  PCa.  Quantitative 
real  time  polymerase  chain  reaction  (QRT-PCR)  was  performed  using  the  7300  Real- 


Time  PCR  System  (Applied  Biosystems,  Foster  City,  CA)  to  validate  and  quantify 
alternative  splicing  events.  Primers  were  designed  to  amplify  the  flanking  regions  of 
skipped  exons  or  the  junctions  across  catenated  exons  of  variant  mRNAs  (see 
illustration  in  Fig,  3b).  Amplified  RT-PCR  products  were  quantified  and  normalized 
to  house-keeping  genes,  EIFIAX and  PPAl,  as  previously  described  Primer 
sequences  for  RT-PCR  validation  are  listed  in  Supplementary  Table  SI. 

Molecular  cloning  of  PIK3CD-S  and  PIK3CD-L  variants  and  overexpression  of 
variants  in  PCa  cell  lines,  RT-PCR  was  performed  to  amplify  PIK3CD-L  and 
PIK3CD-S  transcript  variants  from  purified  RNA  of  PC-3,  VCaP  and  MDA  PCa  2b 
cells  (ATCC,  Manassas,  VA).  PC-3  and  VCaP  cells  were  maintained  DMEM  medium 
(Life  Technologies,  Gaithersburg,  MD)  supplemented  with  10%  fetal  bovine  serum 
(LBS),  while  MDA  PCa  2b  cells  were  grown  in  BRFF-HPC-1  medium  (AthenES, 
Baltimore,  MD)  supplemented  with  20%  FBS.  All  the  cell  lines  were  grown  at  37°C 
and  5%  CO2.  Primers  were  designed  according  to  NCBI  reference  sequences  of 
PIK3CD  mRNA  (NM_005026.3).  The  forward  primer  contained  the  start  codon 
(bold)  (5’-ATGCCCCCTGGGGTGGACT-3’)  and  the  reverse  primer  was  upstream 
of  the  stop  codon  (5’-CTGCCTGTTGTCTTTGGACA-3’).  Pull  length  PCR  products 
were  ligated  into  pcDNA3.1A^5-His  TOPO  vector  (K4800-01,  Invitrogen,  Grand 
Island,  NY)  using  the  manufacturer’s  protocol.  8-10  independent  clones  were  selected 
for  each  of  the  amplified  PIK3CD-L  and  PIK3CD-S  variants  and  sequence  verified. 
The  consensus  sequences  of  PIK3CD-S  and  PIK3CD-L  mRNAs  were  deposited  to 
GeneBank  (Accession  number  KU6121 16  and  KU6121 17).  The  plasmids  pcDNA3.1- 
PIK3CD-LA^5-His  and  pcDNA3.1-PIK3CD-SA^5-Flis  were  individually  transfected 
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into  the  PCa  cell  lines  (VCaP  and  PC-3)  using  the  cationic  lipid  mediated  method  to 
establish  stable  cell  lines  overexpressing  PIK3CD-L  or  PIK3CD-S. 

Functional  analysis  of  PCa  cell  lines  following  transfection  of  slRNAs  targeting 
the  PIK3CD-S  or  the  PIK3CD-L  variant.  VCaP  and  MDA  PCa  2b  cells  were  grown 
in  DMEM  medium  with  10%  FBS  for  24h,  then  were  transfected  for  24h  with 
siRNAs  (50nM)  designed  to  target  splice  variants  of  PIK3CD,  FGFR3,  TSC2  or 
RASGRP2  using  DharmaFECT4  transfection  reagent  (Dharmacon,  Fafayette,  CO), 
according  to  the  manufacturer's  protocol.  In-vitro  functional  assays,  including  cell 
proliferation  and  invasion  were  performed  following  siRNA  transfections  for  24  h. 
Cell  proliferation  and  invasion  assays  were  performed  using  BrdU  Cell  Proliferation 
Assay  kit  (Calbiochem,  Billerica,  MA)  and  the  Matrigel  Invasion  Chambers  (BD 
Biosciences,  San  Jose,  CA),  respectively,  as  previously  described 

Antibodies.  Antibodies  used  in  western  blot  analysis  were  rabbit  monoclonal 
antibodies  for  pAKT^^'^^°^  pAKT^‘=''^’^  AKT,  pmTOR,  mTOR,  pS6  and  S6  (#2965, 
#4058,  #4691,  #2971,  #2983,  #4857  and  #2983,  Cell  Signaling  Technology,  Danvers, 
MA),  rabbit  polyclonal  antibodies  for  His-tag  (abl8184.  Abeam,  Cambridge,  MA), 
PI3K5  and  P-actin  (sc-55589  and  sc-4778,  Santa  Cruz  Biotechnology,  Santa  Cruz, 
CA).  Horseradish  peroxidase  (HRP)-conjugated  secondary  antibodies  for  rabbit  and 
mouse  IgG  were  purchased  from  Southern  Biotech  (Birmingham,  AF). 

In  vivo  xenograft  and  metastasis  models.  All  animal  work  was  approved  by  the 
George  Washington  University  Institutional  Animal  Care  and  Use  Committee 
(protocol  A272).  Four  to  six  week  old,  male  NOD-SCID  mice  were  purchased  from 
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the  Jackson  Laboratory  (Bar  Harbor,  ME).  To  establish  a  PCa  xenograft  model,  2  x 
10^  PC-3  cells  stably  overexpressing  PIK3CD-L  or  PIK3CD-S  were  subcutaneously 
injected  into  the  left  flank  of  NOD-SCID  mice.  Tumor  xenograft  growth  was 
measured  with  calipers  and  the  volume  was  determined  as  1/2  x  length  x  width  . 

Mice  were  randomized  into  groups  once  the  average  tumor  size  reached  -'200  mm^ 
and  treated  with  vehicle  (PBS)  or  CAL-101  (50mg/kg)  through  daily  i.p.  injections. 
After  30  days,  mice  were  sacrificed  and  the  dissected  xenografts  were  photographed 
and  weighed. 

To  establish  the  PCa  metastasis  model,  1x10^  PC-3  cells  stably 
overexpressing  PIK3CD-L  or  PIK3CD-S  were  injected  into  the  tail  vein  of  NOD- 
SCID  mice.  The  mice  were  then  treated  with  vehicle  or  CAL-101  (50mg/kg)  via  i.p. 
injections,  3  times  a  week.  After  8  weeks,  lungs  of  mice  were  harvested  and  stained 
with  India  ink  and  bleached  with  Fekete’s  solution  (70%  ethanol,  3.7%  formaldehyde, 
0.75  M  glacial  acetic  acid).  India  ink-stained  lungs  were  photographed  and  lung 
metastases  were  quantified  using  the  NIH  Image!  program 

Purification  of  His-tagged  PI3K5  protein.  PC-3  cells  stably  overexpressing 
PIK3CD-L  or  PIK3CD-S  were  maintained  in  DMEM  media  (Life  Technologies, 
Gaithersburg,  MD)  supplemented  with  10%  fetal  bovine  serum  (FBS).  After  growing 
the  cells  for  24  h,  cell  extracts  were  prepared  and  His-tagged  PI3K5  protein  was 
purified  using  a  column  HisPur  Ni-NTA  purification  kit  (Pierce  Biotechnology, 
Rockford,  IE).  Briefly,  cell  lysates  were  mixed  with  Ni-NTA  resin  and  incubated  at 
room  temperature  for  30  min.  After  incubation,  the  resin  was  washed  with  wash 
buffer  (25mM  imidazole,  pH  7.4)  and  applied  to  a  HisPur  Ni-NTA  spin  column, 
centrifuged  and  wash  buffer  eluate  discarded  after  centrifugations.  His-tagged 
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proteins  were  eluted  from  the  resin  by  adding  one-resin-bed  volume  of  elution  buffer 
(250mM  imidazole,  pH  7.4).  The  purified  PI3K6-His  protein  was  mixed  with  2x 
Laemmli  sample  buffer,  boiled  and  analyzed  by  immunoblotting. 

Co-immunoprecipitation  (Co-IP)  of  PI3K6/p85  complex.  Plasmids  pcDNAS.l- 
PIKSCD-Sm-His  (or  poDNA3.1-PIK3CD-LA^5-His)  andpSV-p85a  (Addgene, 
Cambridge,  MA)  were  eo-transfeeted  into  PC-3  cells.  After  growing  the  cells  for  48 
h,  the  co-transfected  cells  were  harvested  and  cells  were  lysed  with  RIPA  lysis  buffer 
(Santa  Cruz  Biotechnology,  Santa  Cruz,  CA).  The  cell  lysates  were  then  subjected  to 
Co-IP  assays  with  anti-His  antibody  (abl8184.  Abeam,  Cambridge,  MA)  and 
immobilized  on  protein  G-Sepharose  beads  (Thermo  Scientific,  Waltham,  MA).  Cell 
lysates  and  precipitates  were  subjected  to  western  blotting,  and  visualized  by 
enhanced  chemiluminescence  (ECL)  system  (Thermo  Scientific,  Waltham,  MA). 

In  vitro  assay  of  PI3K5  activity,  PI3K5  activity  was  evaluated  with  a  PI3K 
activity/inhibitor  assay  kit  (Millipore,  Billerica,  MA)  according  to  the  manufacturer’s 
instructions.  Briefly,  purified  His-tagged  PI3K5-L  or  PI3K5-S  isoform  was  pretreated 
with  the  PI3K6  inhibitor  (lOOnM  of  wortmannin  or  lOOnM  of  CAL- 101)  or  vehicle  in 
96  well  plates  for  10  min  and  subjected  to  a  competitive  ELISA.  PIP2  substrate  and 
kinase  reaction  buffer  were  added  to  the  pretreated  His-tagged  PI3K6-L  or  PI3K5-S 
isoform  and  incubated  at  room  temperature  for  1  h.  After  incubation,  biotinylated 
PIP3  and  GST-GRPl  working  solutions  were  added  to  the  wells  and  the  reaction 
samples  were  further  incubated  at  room  temperature  for  1  h.  Plates  were  washed  three 
times  with  IxTBST  (50  mM  Tris-Cl,  pH  7.5.  150  mM  NaCl)  and  incubated  with  SA- 
HRP  (Streptavidin-Horseradish  Peroxidase  conjugate,  1 .25  mg/ml)  at  room 
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temperature  for  1  h.  After  incubation,  plates  were  washed  three  times  and  incubated 
with  100  pi  of  TMB  (3,3’,5,5’-tetramethylbenzidine,  1  mg/ml)  substrate  solution  at 
room  temperature  for  5-20  min.  Reactions  were  stopped  by  adding  100  pi  of  stop 
solution  and  plates  read  at  450nm.  The  colorimetric  signal  was  inversely  proportional 
to  the  amount  of  PIP3  produced  by  PI3K  activity  and  the  relative  amount  of  PIP3 
produced  was  determined  with  a  standard  curve. 

Additional  methods  can  be  accessed  in  the  Supplementary  Methods  section. 
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Fig,  1,  Heatmaps  of  alternative  splicing  differences  in  AA  PCa  vs,  EA  PCa  and 

AA  NP  vs,  EA  NP  specimens,  (a)  PCA  plot  and  2D-clustergram  depicting  3,112 
significant  differentially  expressed  exons  in  20  independent  AA  PCa  vs.  15 
independent  EA  PCa  specimens,  (b)  PCA  plot  and  2D-clustergram  depicting  3,384 
significant  differentially  expressed  exons  in  20  AA  NP  vs.  15  EA  NP  specimens.  AA 
and  EA  specimens  are  represented  by  red  and  blue  circles/bars,  respectively.  Rows 
represent  specimens  and  columns  represent  exons  in  hierarchical  clustergrams.  Eog2 
expression  values  of  exons  were  subjected  to  2D  unsupervised  hierarchical  clustering 
using  average  linkage  method  and  Euclidean  distance,  (c)  Venn  diagrams  of  DS 
events  in  AA  PCa  vs.  EA  PCa  and  AA  NP  vs.  EA  NP,  DS  events  in  EA  PCa  vs.  EA 
NP  and  AA  PCa  vs.  AA  NP  and  differentially  expressed  genes  in  AA  PCa  vs.  EA 
PCa  and  DS  events  in  AA  PCa  vs.  EA  PCa.  d)  Majority  of  the  genes  with  DS  events 
in  AA  PCa  vs.  EA  PCa  were  functionally  associated  with  cancer.  Top  3  ‘other 
diseases’  were  gastrointestinal  disease,  organismal  injury  and  abnormalities  and 
reproductive  system  disease. 
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Fig.  2.  Representation  of  differential  splicing  events  in  a  composite  oncogenic 
signaling  pathway,  DS  events  in  AA  PCa  vs.  EA  PCa  were  frequently  deteeted  in 
multiple  oncogenic  signaling  pathways.  Depicted  pathway  is  a  composite  of 
PBK/AKT/mTOR,  RAS/RAF/MAPK,  CDK/RBl,  MDM2/TP53  and 
WNT/GSK3/CTNNB1/APC  signaling.  Open  figures  indicate  oncogenes  and  closed 
figures  indicate  tumor  suppressor  genes.  “DS”  indicates  that  a  differential  alternative 
splicing  event  was  detected  for  a  particular  oncogene  or  tumor  suppressor  gene. 
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Fig,  3,  Validation  of  differential  alternative  splicing  of  oncogenes  and  tumor 
suppressor  genes  in  AA  vs  EA  PCa  specimens,  (a)  Schematic  representation  of  DS 
events  within  the  indicated  genes  in  AA  PCa  (top  dashed  lines  connecting  exons)  and 
EA  PCa  (bottom  dashed  lines)  based  on  Alternative  Splicing  ANOVA  model.  Closed 
arrowheads  below  exons  represent  primer  location  in  QRT-PCR  validation  of 
alternatively  spliced  transcripts,  (b)  Representative  RT-PCR  results  validating  race- 
specificZ-enriched  variant  transcripts  in  either  AA  PCa  or  EA  PCa  specimens.  Shown 
are  the  RT-PCR  results  for  the  AA-specific/-enriched  variants  PI3KCD-S,  FGFR3-S, 
TSC2-S,  ITGA4-L,  MET-L,  NFl-L,  BAKl-L  and  RASGRP2-b;  and  EA-specificZ- 
enriched  variants  PIK3CD-L,  FGFR3-L,  TSC2-L,  ITGA4-S,  MET-S,  NFl-S,  BAKl-S 
and  RASGRP2-a.  Each  lane  represents  an  RT-PCR  result  from  an  independent  PCa 
specimen  that  was  also  interrogated  in  exon  array  experiments.  RT-PCR  of  EIFIAX 
and  PPAl  transcripts  served  as  loading  controls.  QRT-PCR  results  is  summarized  in 
Supplementary  Fig,  S4.  (c)  Additional  validation  and  quantification  of  the  ratio  of 
the  AA-enriched  PI3KCD-S  (short)  variant  and  the  race-independent  PI3KCD-L 
(long)  variant  in  a  separate  cohort  of  PCa  patient  specimens.  RNA  was  isolated  from 
n=32  AA  PCa  and  n=30  EA  PCa  specimens  and  subjected  to  QRT-PCR.  Shown  is  a 
plot  of  the  ratio  of  S/E.  EIFIAX  and  PPAl  transcripts  served  as  internal  normalization 
controls.  *  E  <  0.05  using  Student  t-test. 
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Fig.  4.  Selective  knockdown  of  the  AA-enriched  PIK3CD-S  variant  or  race- 
independent  PIK3CD-L  variant  has  opposing  effects  on  proliferation,  invasion 
and  AKT/mTOR  signaling  in  AA  PCa  cells,  (a)  RT-PCR  of  PIK3CD-L  and 
PIK3CD-S  variants  following  knockdown  in  the  EA  PCa  cell  line  VCaP  and  AA  PCa 
cell  line  MDA  PCa  2b.  Specific  knockdown  of  the  PIK3CD-L  variant  (left  panel)  was 
accomplished  with  an  exon  20-targeting  siRNA  (siP2o),  while  knockdown  of  the 
PIK3CD-S  variant  (right  panel)  was  achieved  with  an  siRNA  targeting  the  region 
spanning  exons  19  and  21  (siPj).  Closed  arrowheads  below  exons  represent  primer 
location  for  QRT-PCR  validation  of  alternatively  spliced  transcript  knockdowns. 
Knockdown  efficiency  of  siP2o  and  siPj  siRNAs  was  determined  by  the  S/L  ratio 
derived  from  RT-PCR  reactions,  ns  =  nonsense  siRNA  treatment.  Representative 
images  of  n  =  4  independent  knockdown  experiments,  (b)  Proliferation  and  invasion 
of  VCaP  and  MDA  PCa  2b  following  knockdown  of  the  PIK3CD-S  (left  panel)  or 
PIK3CD-L  variant  (right  panel).  Data  presented  as  mean  ±  SEM  of  n=3-5  independent 
experiments  for  each  treatment  group.  *7*  <  0.05  by  ANOVA  with  post-hoc  Tukey. 

(c)  Western  blot  analysis  of  AKT/mTOR  signaling  following  knockdown  of  the 
PIK3CD-L  (top  panel)  or  PIK3CD-S  variant  (bottom  panel)  in  VCaP  and  MDA  PCa 
2b.  Eevel  of  AKT,  mTOR  and  ribosomal  S6  kinase  activities  is  reflected  by  the 
amount  of  phospho-AKT  (pAKT),  phospho-mTOR  (pmTOR)  and  phospho-ribosomal 
protein  S6  (pS6)  immunoblotting,  respectively.  (3-actin  served  as  loading  control. 
Representative  images  from  n  =  3-4  independent  western  blot  experiments. 
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Fig,  5,  PI3K6-S,  but  not  PI3K6-L,  is  resistant  to  small  molecule  inhibition  of 
PI3K/AKT/mTOR  signaling  and  proliferation,  (a)  Assessment  of 
PBK/AKT/mTOR  signaling  following  treatment  with  vehicle  (saline)  or  CAL-101 
(100  nM,  24  hrs).  PBK/AKT/mTOR  signaling  was  assessed  by  western  blot  analysis 
with  phospho-antibodies  to  AKT  (pAKT)  and  S6  ribosomal  protein  (pS6).  (3-actin 
served  as  a  loading  control.  His-tag  antibody  was  used  to  demonstrate  equal 
expression  of  His-tagged  variant  PBK5  protein  in  stably  transfected  cell  lines. 
Representative  images  from  n  =  3-4  independent  western  blot  experiments. 

*  Significantly  different  from  corresponding  vehicle  control,  P  <  0.05  by  ANOVA 
with  Dunnett’s  post-hoc  test,  (b)  Proliferation  in  VCaP  and  PC-3  cells  stably 
overexpressing  the  PIK3CD-S  variant  or  PIK3CD-L  variant  following  treatment  with 
vehicle  (saline)  or  selective  PBK5  small  molecule  inhibitor  CAL-101  (24  hrs),  and 
(c)  treatment  with  vehicle  (saline)  or  selective  AKT  small  molecule  inhibitor  MK- 
2206  (24  hrs).  Proliferation  was  assessed  using  a  BrdU  labeling  assay.  Data  presented 
as  the  mean  ±  SEM  of  n=4-6  independent  experiments  for  each  treatment  group.  *  or 
#Significantly  different  from  corresponding  vehicle  control,  P  <  0.05  by  ANOVA 
with  Dunnett’s  post-hoc  test. 
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Fig,  6,  PI3K6-S,  but  not  PI3K6-L,  is  resistant  to  small  molecule  inhibition  of 
xenograft  growth  and  metastasis,  (a)  Growth  of  PC-3  cells  stably  overexpressing 
the  PIK3CD-S  variant  in  NOD-SCID  mice  is  resistant  to  CAL-101  treatment 
(50mg/kg  i.p.  3  times  a  week).  In  eontrast,  growth  of  PC-3  cells  stably  overexpressing 
the  PIK3CD-L  variant  is  sensitive  to  CAL-101  treatment.  Data  represent  the  mean 
tumor  size  ±  SEM  of  n=10  independent  miee  for  each  treatment  group  at  eaeh  time 
point.  *Signiiioantly  different  from  saline  treated  group,  P  <  0.05  by  ANOVA  with 
post-hoc  Tukey.  (b)  Tumor  weights  and  gross  morphology  of  tumor  xenografts  from 
(a).  Box- and- whisker  plot  represents  mean  xenograft  weight  in  miee  after  30-day 
vehiele  or  CAL-101  treatment.  *P  <  0.05  by  ANOVA  with  Dunnett’s  post-hoe  test; 
n=10  independent  mice  for  each  treatment  group,  (e)  Quantifieation  of  lung 
metastases  in  NOD-SCID  miee.  PC-3  cells  stably  overexpressing  PIK3CD-L  or 
PIK3CD-S  were  injeeted  into  the  tail  vein  of  NOD-SCID  mice  treated  with  vehiele  or 
CAL-101  (50mg/kg  i.p.  3  times  a  week).  After  8  weeks,  lungs  were  harvested  and 
stained  with  India  ink  and  bleached  with  Fekete’s  solution  for  visualization  of 
metastatic  nodules  (white-eolored  areas).  Data  presented  as  mean  ±  SEM  of  n  =  10  for 
each  treatment  group.  *P  <  0.05  by  ANOVA  with  Dunnett’s  post-hoc  test,  (d) 
Representative  India  ink-stained  lungs  from  treatment  groups  analyzed  in  (e). 
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Fig,  7,  Cell-free  kinase  assay  of  PI3K6  isoforms  and  tertiary  structure  modeling 
of  PIK36-small  molecule  inhibitor  interactions,  (a)  Schematic  representation  of 
protein  domains  of  PI3K5-L  and  -S  isoforms.  Adaptor  (p85)-binding,  RAS  binding, 
C2,  helix  and  catalytic  domains  are  highlighted.  Phosphorylation  and  ubiquitination 
sites  (S,  T,  Y  and  K)  and  the  region  encoded  by  exon  20  (56  amino  acids)  residing  in 
the  eatalytic  domain  are  indieated.  (b)  Co-immunoprecipitation  (Co-IP)  of  PI3K6/p85 
complex  followed  by  western  blot,  (c)  and  PI3K  activity  assays.  S,  supernatant;  W, 
wash  fraction;  E,  eluted  fraction.  Blots  are  representative  of  3-5  independent 
experiments  with  similar  results.  Anti-His  antibody  was  used  in  the  Co-IP 
experiments  and  anti-His,  anti-p85a  and  anti-aetin  antibodies  were  used  in  the 
western  blotting.  *  Significantly  different  kinase  activities  in  total  lysates  of  PI3K5-S 
vs.  PI3K5-L  expressing  cells,  or  purified  PI3K5-S/p85  vs.  PI3K5-L/p85  complexes.  P 
<0.05  using  Student  t-test.  Data  presented  as  mean  ±  SEM  of  n  =  4  for  eaeh  treatment 
group,  (d)  Purification  of  His-tagged  PI3K5-E  and  -S  isoforms.  Western  blot  analysis 
of  Ni-NTA  resin-purified  PI3K5  isoforms  from  transfected  PC-3  cells  using  His  and 
PI3K5  antibodies.  FT,  flow-through;  W,  wash  fraction;  E,  eluted  fraction.  Blots  are 
representative  of  4  independent  experiments  with  similar  results.  Closed  arrowheads 
indieate  PI3K5  isoforms,  (e)  Cell-free  kinase  assay  of  L  and  S  isoforms  of  PI3K5  in 
the  presence  of  vehicle  (PBS),  lOOnM  wortmannin  or  lOOnM  CAL-101. 

*Significantly  different  from  vehicle  control-treated  PI3K5-L  isoform.  P  <  0.05  by 
ANOVA  with  Dunnetf  s  post-hoc  test.  Data  presented  as  mean  ±  SEM  of  n  =  4-5  for 
each  treatment  group. 
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